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ABSTRACT 

Various methods have been introduced for dynamic braking 
of ships. Depending on the prime mover and transmission 
system, dynamic braking is accomplished by electric power 
dissipation (in electric driven ships) or clutches and brakes, 
or air compression (in diesel engine propulsion systems) . 

In this study a gas turbine - electric drive combination 
is used as the means of propulsion; the usual method for 
dynamic braking is the electric power dissipation method 
using resistors. 

This thesis is concerned with an investigation of dif- 
ferent methods for dynamic braking in an attempt to find a 
way to decrease the size of the braking resistors needed and 
their associated equipment, or to eliminate the need for these 
by using another type of dynamic braking. 
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I. 



INTRODUCTION 



This thesis is concerned with the study of a ship's 
dynamic braking and reversing system using a combination of 
gas turbine - electric propulsion motors. 

Because of the ever present danger of collisions, the 
problem of crash-stopping a ship is becoming more urgent both 
in the Navy and the Merchant Marine. 

The use of a reduction gear mechanism for gas turbine 
combatant ships is quite attractive but in this case, during 
the reversing phase it introduces a complex reduction gear 
mechanism requiring large clutches for the accomplishment of 
forward and astern motion; thus, an electric drive is also 
attractive. In this study the electric drive has been assumed 
and various methods for dynamic braking are investigated. 

The first part of the thesis has been devoted to the 
study of different open loop systems in order to select the 
one that best suits the required specifications; these have 
been defined as follows: 

a. Stop the ship in minimum time. 

b. Try to minimize additional equipment for dynamic 
braking. 

The second part of the thesis is concerned with the 
study of the control system selected in the first part imple- 
mented by a microprocessor-based controller. 
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II. THE MODEL SHIP 



For the simulation of the problem in the computer (IBM 
360), a model ship was derived. For this, a destroyer escort 
type vessel was selected with an approximate displacement of 
4000 tons, a length of 400 feet and a top speed of 30 knots 
[Ref. 1] . 

A single fixed propeller was selected and the propulsion 
plant consists of a marine gas turbine (FT4A-2) , which directly 
drives the main propulsion generator to produce the necessary 
electric power for the main motor which drives the propeller 
shaft. A schematic diagram of the propulsion system is shown 
in Figure 1. 

The simulation of the ship and propulsion plant uses the 
above model with somewhat idealized hydrodynamic and propulsion 
plant dynamics. 

Since the equations governing the response of the ship 
and propulsion plant are complicated non-linear differential 
equations, the digital computer was used to solve the equa- 
tions using numerical analysis methods. 

As a consequence of the idealized model the following 
non-linear functions were stored in the computer memory as 
look up tables during the simulation: 

a. Thrust reduction coefficient 1 - t 

b. Wake fraction coefficient 1 - W 

c. Ship's resistance R 
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d. Propeller thrust coefficient 
and e. Propeller torque coefficient 

where the thrust reduction coefficient (1-t) is a measure of 
the propeller and hull interaction, and the wake fraction co- 
efficient (1-W) is a measure of the water velocity near the 
propeller; the propeller thrust (C_) , and torque (C ) coef- 
ficients are measures of the efficiency of the propeller. 

All the above parameters, along with the corresponding 
computer programs, are given in Appendix B. 
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Fig.l- Schematic representation of the propulsion 
system 
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III. DYNAMIC BRAKING 



The propeller speed can be reduced quickly from some 
initial large valve during the coast down phase by dynamic 
braking. 

Depending on the prime mover and the transmission system, 
dynamic braking is accomplished by electric power dissipation 
(with electric drive) , clutches (reversing gear) , air com- 
pression (diesel engine) , or with a water or air brake on the 
shaft. 

During dynamic braking, the power delivered by the drive 
train must be equal to the power absorbed by the braking system 
and the losses. Since the drive train with its inertia is 
being decelerated rapidly, there is a large "transient inertial 
horsepower" due to the shaft and propeller speed deceleration, 
as well as rotational power. 

In terms of a torque equation, this can be expressed as: 



2 tt t + = 



dt 



Q mot " Q p Q fr ” °BR 



( 1 ) 



where, Q. 



mot 



FR 



brake 



is the torque developed by the electric motor on 
the propeller shaft 

is the torque developed on the propeller by the 
mass of water passing through its blades due to 
the motion of the ship. 

is the friction torque on the shaft and motor 
is the torque developed by the dynamic braking 
device, opposing the motion of the propeller shaft 
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is the total moment of inertia of the motor rotor. 



propeller shaft and propeller (i.e., I 



N 



I + I ) 

shf t p 

is the speed of the propeller. 



I + 
m 



For this study the following means of dynamic braking 
will be considered: 

a. Electric power dissipation 

b. An air brake 

c. Combinations of the above 

d. Use of the inertia of moving parts as load and 
absorbing medium. 

All runs have assumed an initial maximum speed of 30 
knots and a maximum propeller speed of 230 rpm. 
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IV. SHIP'S MOTION EQUATIONS 



A moving ship is a body with six degrees of freedom. 
These degrees of freedom are generally chosen as follows: 

a. Linear displacement along the three axes through the 
center of gravity. 

b. Rotations around the three axes through the center 
of gravity. 

If we assume ship's motion in calm water with no turning 
maneuvers, then the motion of the ship is best described by 
Newton's law of motion for the ship in translation and the 
propeller in rotation [2], as follows: 

thrust equation 

M ~=(l-t).z .T (V ,n)-R(V) (2) 

dt P P P 



and torque equation 



<w 



- «fr -<5, 



(3) 



where M is the mass of the ship 

V is the velocity of the ship 

Z is the number of shafts 
P 

T (V ) is the propeller thrust as a function of the 
^ ^' n propeller advance speed (V ) and propeller 

speed (N) p 

and R is the ship's resistance as a function of speed, 

The input data required to simulate an acceleration or 
deceleration manuever is developed from the characteristics 
of the ship and propulsion system being modelled, and may be 
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summarized as follows: 



Ship's displacement 
Propeller parameters 

Non-dimentional propeller performance data 
Prime mover torque characteristics 
Polar moment of inertia of rotating system 
The coefficient of added mass (c) , representing the mass of 
water entrained by the ship, has been assumed to be 0.08 in 
all cases as recommended in Ref. 12. 

Ship's resistance (R) has been assumed as a function of 
ship's speed. In addition, the hull-propeller interaction 
is also assumed to be a function of ship speed. The wake 
factor (1-W) and the thrust, reduction factor (1-t) are gener- 
ally available from self-propelled model tests. The relative 
rotational efficiency is assumed to be constant and equal to 
unity. The net thrust applied to the ship is the product of 

the propeller thrust (T ) and the thrust reduction factor 

P 

(1-t) . The ship speed (V) is multiplied by the wake factor 

(1-W) to give the velocity of advance of the propeller (V ) . 

P 

In terms of equations the above are given as follows 

[ 12 ]: 

Propeller thrust equation 



T =(l-t).T 
P a 



( 4 ) 



where T^ is the propeller theoretical thrust 
and 

Velocity of advance of propeller 
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V =(l-w).V 
P 



( 5 ) 



where V is the velocity of the ship. 

The data required for adequate representation of fixed 
pitch propeller performance, throughout the range of operating 
conditions encountered during ship acceleration and deceleration 
maneuvers have been developed by Miniovitch [2] . Miniovitch 
has published dimentionless performance data for three-bladed 
propellers with blade area ratios from 0.6 to 1.6 in incre- 
ments of 0.2. 

The data are presented in the form of thrust (k^) and 
torque (k ) coefficients versus the advance coefficient (J) 

q 

2 2 

and also in the form of K /J and K /J vs 1/J. The modified 

t q 

2 2 

thrust (k /J ) and torque (k /J ) coefficients versus the 

U SI 

reciprocal of the advance coefficient (1/J) are required in 
order to account for the case where the propeller rotational 
speed is very low or equal to zero. The data extends over 
the four distinct operating conditions to be encountered by 
the propeller, namely: 

Ship travelling ahead, propeller turning ahead 
Ship travelling ahead, propeller turning astern 
Ship travelling astern, propeller turning astern 
Ship travelling astern, propeller turning ahead. 

The data for a three bladed propeller have been stored in 
the computer as a look up table . 
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The above, in terms of equations can be expressed as 



follows : 



Propeller advance coefficient 



o ?' 1/2 

J=N .D/(V + (N .D) ) 



P* P P 



( 6 ) 



Propeller thrust 



T =C rr .p.D 2 . (V 2 +(N .D) 2 ) 
a r ^ p p 



(7) 



and Propeller torque 



Q =C .p.D 3 . (V 2 +(N .D) 2 ) 

p q p p 



( 8 ) 



where 

C T is the propeller thrust coefficient 
Cg is the propeller torque coefficient 
D is the propeller diameter 
p is the sea water density 
Np is the propeller speed 

The polar moment of inertia of the rotating system has been 
calculated in Appendix A. The inertia of the water mass 
entrained by the propeller has also been calculated. 

If there is more than one propulsion shaft, then Eq. (3) 
will be repeated for each additional shaft [2] . 

The direction of rotation of the free turbine is always 
clockwise as is the direction of rotation of the generator 
set. The direction of rotation of the motor and the propeller 
shaft is assumed positive if the shaft rotates clockwise when 
viewed from the shaft end. 
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Figure 2 shows the assumed rotation and the relation of 
the different torques acting on the propeller shaft for for- 
ward motion of the ship. 



The frictional torque of the propeller shaft (Qp R ) has been 
assumed constant throughout the simulation; actually it is a 
function of the shaft speed and can be represented also as a 
stored table look up [7] . The frictional torques of the 
motor and generator and their windage torques have been 
assumed negligible. 

✓ 

A. THE GAS TURBINE-GENERATOR SYSTEM 

The gas turbine-generator set speed is controlled by the 
influence of the engine torque, the generator load torque and 





Fig. 2- Assumed notation on application of torques 
and direction of rotation. 





(9) 
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the friction losses of the system. This can be expressed 
as follows [6] : 



2r».I 



.^t= Q 

tot dt° V E 




(10) 



where 



Q_ is the gas turbine drive torque. It is a non- 
E 

linear function of fuel consumption (W_) and 

r 

gas turbine speed [8], 

Qgen is the generator load torque given by [4]: 



Q = K .i p (t) .i (t) 
y gen g fg a 



( 11 ) 



where 

Kg' is the torque constant of the generator 

i^g(t) is the generator field current 

i (t) is the armature current 
a 

I. . is the sum of inertias of motor rotor 
tot 

(I ) and gas turbines power turbine (I . ) 
gen gt 

Q^ r is the friction torque which is made up of 
brush friction and windage components, given 
by [4] : 



where 



Q- = Q, + k .Nj: 
x fr x br w gt 

Q^ r is the brush torque 

is the windage coefficient 



( 12 ) 



The engine is set to run as a nominally constant speed engine 
with a proportional governor. Since the output torque of the 
gas turbine is a non-linear function of fuel consumption and 
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gas-turbine speed, data for the output torque have also been 
stored in the computer as a look up table. 



B. ELECTRICAL TRANSMISSION SYSTEM 

The back emf generated by the generator and the motor 
are given by [4] : 



E = K .i p (t) .N 
g g fg gt 



(13) 



and 



E =K .ip(t) .N 
m m rm p 



(14) 



where 

i,. , i,. are the field currents of motor and gen- 
fm fg 

erator, respectively 
Kg is the generator electrical constant 
Km is the motor electrical constant 
Em is the back emf of the motor 
Eg is the back emf of the generator 
Ignoring the effects of the armature inductances of motor and 
generator, for maneuvering conditions the armature current 
will be given by 

i (t)=(E (t)-E (t) )/R. (15) 

a a m t 



where 




R + R 



m g 



(16) 



and where 

Rm is the motor armature resistance 
Rg is the generator armature resistance. 
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V. DYNAMIC BRAKING USING ELECTRIC POWER DISSIPATION METHODS 



A braking resistor is utilized here in order to dissipate 
the power produced by the motor due to the motion of the shaft. 

When the CRASH-STOP command is given, the main generator 
is disconnected from the system and the motor is left to rotate 
due to its inertia plus the applied torque on the propeller 
due to the moving mass of water through its blades due to the 
motion of the ship. 

Thus, the motor effectively becomes a generator; by 
leaving its terminals open, since no current is circulating, 
little opposing torque will be produced. If, on the other 
hand, a resistor is connected at the terminals of the motor, 
a current will circulate and thus a torque opposing the motion 
of the motor will be produced and the propeller speed will 
drop. 

In order to have maximum opposing torque on the propeller 
shaft, the maintenance of maximum circulating armature current 
is required; but since the voltage generated at the terminal 
of the motor - when functioning as a generator - is variable, 
the above condition can only be met by a variable braking 
resistor or by a variable motor field current. 

In this study only the variable resistor method has been 
investigated. Figure 3 shows a schematic representation of 
the propulsion system after the CRASH-STOP command has been 
initiated. 



27 



The gas turbine can be held at load, if desired, with 

the simultaneous addition of a fixed resistor R as shown in 

o 

Figure 3. 

The equivalent circuit of the motor is shown in Figure 
4 after t = 0, i.e., after the CRASH-STOP command has been 
issued. 

Applying Kirchoff's law around the loop we get 



V (t)= E (t)-i (t) .R -L ^a (t) 
m m a m m dt 



(17) 



where 

Vm is the voltage generated at the motor terminals 
and Lm is the motor armature inductance 
Due to the addition of the braking resistor the motor now is 
• effectively working as a generator; the load being the braking 
resistor Rg and the prime mover being the torque produced 
on the propeller due to the motion of sea water through its 
blades . 

We already have assumed that the reactance of the arma- 
ture circuit is negligible so that Eq. (17) is transformed to 



V =E (t)-i (t).R m (18) 

mm a. m 

Applying Newton's law for the motor-propeller system, we get 
Eq. (19) , i.e. , 

2TT,(I mot +I L +I shft ) "dtP =( ^iiot" '^fr (19 

where I . = moment of inertia of motor 

mot 

I T = moment of inertia of propeller 
J-i 

I shft = moraent inertia of propulsion shaft 
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Fig. 3. Schematic representation of the 

propulsion system after disconnection (t 0) 




Fig. 4. Equivalent circuit representation of 

the separately excited DC motor at t 0 
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xUJo 



and 



the torque produced by the motor is given by 



Q 



mot ' 






= K' . ij- (t) .i (t) 
m rm a 



( 20 ) 



where K ' is the torque constant of the motor, 

m 

As can be seen from Eq. (20), since the armature current 
has opposite direction, Q mot will be negative as far as i fm (t) 
remains positive, i.e., as long as we do not reverse the motor. 
Therefore, the right hand side of Eq. (19) is negative and 
the motor decelerates faster than in the case where the ter- 
minals were kept open (i.e., Q mot = 0 ) • 

From Eq. (20) it can also be seen that the bigger Q mot 
is, the faster the deceleration of the propeller shaft will 
be. Thus the armature current must be kept at maximum rated 
value which can be accomplished by varying the braking resistor 



V 

In order to calculate the rate of change of the braking 
resistor Rg, for maximum deceleration, the program of Appendix 
B was used with the armature current held constant at its 
maximum rated value. The values of the braking resistor can 
then be calculated from Eq. (18) where solving for Rg we get 



p _JL(t)-i .R m 

Rg m — — — 3. — ni- 



( 21 ) 



where 



v =R d .I 
m B c 



( 22 ) 



The calculated values of the braking resistor are shown 
in Table I and have been plotted versus time in Figure 11. 
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As can be seen the law of variation of the braking resistor 
is not linear, therefore an approximation was made by fitting 
straight lines to the calculated values. 

The values of the linearized braking resistor are given 
in Table I also, and they are plotted versus time in Figure 
12 . 

At the instant the shaft speed reaches a certain value"*" 
determined by the maximum loading capabilities of the system 
and assuming that the generator field current has been reversed, 
the system is reconnected as before and astern maneuvers are 
applied. 

The armature current in the loop will be given by 

i «st (t)=tE g (t> -y t))/( v R g ) <23) 

where Eg(t) is the back emf of the generator given by Eq. (13) 
and Em(t) is the back emf of the motor given by Eq. (14) . Rm 
and Rg are the armature resistances of motor and generator 
respectively . 

The results of this method of dynamic braking using a 
linearized braking resistor are shown in Figure 13 through 
Figure 16, along with results that would have been obtained 

if the theoretical variable resistor was used. 

/ 

In particular Figure 13 shows the variation of ship's 
velocity versus time. Due to the fact that the values of the 
linearized braking resistor are quite close with the values 



1 



This value has been calculated in Appendix A to be 13 rpm. 
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of the theoretical one, the differences of speed variation 
are negligible. 

In Figure 14 the variation of propeller speed is shown 
versus time for both cases. The abrupt change of the curve 
slope during dynamic braking, (for 0<i<1.4 sec.) is noticeable 
compared with the change of slope during astern manuevers. 

In Figure 15 and Figure 16 a power dissipation diagram 
is presented. In particular in Figure 15 the power dissipated 
in the braking resistor (PRB) and the internal resistance of 
the motor (PRM) have been plotted against time for the case 
where the theoretical variable braking resistor is used. Also, 
on the same graph, the percentage of dissipated power in the 
braking resistor and the internal resistance have been plotted 
(PERC1 and PERC2 respectively) . 

On the other hand, in Figure 16 the same values but for 
the case of the linearized braking resistor have been plotted. 
The difference in the percentage of dissipated power can be 
noticed since now a linear resistor is being used. Due to 
linearization, less power on a percentage basis is dissipated 
by the braking resistor. Thus, the percentage of power dis- 
sipated on the internal resistance of the motor increases. 

Due to the fact that with a linearized variable braking 
resistor the percentage of power dissipation in the internal 
resistance of the motor increases, another method for dynamic 
braking using a fixed braking resistor was selected. 

In this case, as it can be expected, the time required 
to stop the propeller and eventually the ship,, increases . 

This is because the current circulating in the loop is no 
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longer maximum; therefore the opposing torque generated by 
the motor (given by Eq. 20) , decreases along with the current. 

The time required to stop the propeller with a variable 
resistor was found to be 1.8 seconds while with a fixed re- 
sistor this time increases to 28 seconds. 

The results of the method where we use a fixed braking 
resistor are shown in Figure 17 to Figure 19. In particular 
Figure 17 shows the variation of ship speed versus time and 
it can be noticed that the ship's speed at t = 30 sec is still 
above 12 knots while from Figure 14 the time to stop the ship 
when a variable resistor is used, was found to be 23.2 seconds. 

In Figure 18 the variation of propeller speed versus time 
has been plotted. The largest deceleration rate of the pro- 
peller speed is at the beginning of the dynamic braking phase, 
since then the circulating current is still maximum. But as 
time advances, the circulating current drops, therefore the 
opposing torque produced on the propeller shaft becomes less, 
eventually leaving the system in a feathering propeller con- 
dition. At t = 28 sec the slope changes due to the application 
of the astern motion on the motor. 

In Figure 19 as in Figures 15 and 16, a power dissipation 
diagram is presented. The constant percentage dissipation 
curves PERC1 and PERC2 (for the braking resistor and the 
internal resistance of the motor respectively) can be noticed. 
They indicate that the percentage power dissipated in the 
armature resistor is constant at 3% while the remaining 97% 
is being dissipated in the fixed braking resistor. 
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Table I. Theoretical and Approximated Values of the Braking 
Resistor RB (ohms) 

^ IME THEORETICAL APPROXIMATED 



0.0 


1.72 10 -2 


1.7 10 


-2 
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1.7136 


IQ" 2 


1.6015 


IO -2 
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1.5973 


io -2 


1.5029 


IO -2 
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1.3819 


io -2 


1.4044 


10 -2 
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io -2 


1.3059 


10~ 2 
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io- 2 
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io' 2 
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10 -2 
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io -2 
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The CRASH-STOP command was executed in 23.2 seconds when 
a variable linear braking resistor was used and in t » 35 sec 
when a fixed braking resistor was used. 
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VI. DYNAMIC BRAKING USING MECHANICAL POWER 
DISSIPATION METHODS 



Since the scope of this thesis is to try to minimize 
the additional braking devices required for dynamic braking, 
it may seem strange to consider the introduction of a mechani- 
cal braking device. 

But, on the other hand, if the sizes of the resistors 
found by the application of the previous method are kept in 
mind, associated with the problems encountered with heat dis- 
sipation then it will be found that the use of a mechanical 
brake may help in reducing the size of the resistors. 

There is a large variety of mechanical brakes which can 
be used, such as water brakes, friction brakes, or eddy current 
type brakes. 

Clutches and brakes are used to control the coupling be- 
tween two shafts by means of relatively low power signals. 

This coupling can be either of the continuous type, where the 
degree of coupling varies smoothly with the magnitude of the 
control signal, or of the on-off type where the coupling is 
either zero or maximum. 

A brake is not a torque converter but a transmitter of 
torque. The brake controls speed by permitting slip with 
respect to the prime mover shaft. To accomplish this, the 
brake must absorb energy and thereby generate heat [3] . 

Heating generally establishes the maximum power level 
that can be controlled by the brake by limiting the maximum 
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slip. An analysis of the heating problem must include an 
examination of the duty cycle and the inertial and/or dissi- 
pative characteristics of the load. Heat dissipation is of 
maximum concern in brakes that require the slipping of friction 
shoes where heat is localized and wear is high. 

Since no data were available on a real brake, an air- 
clutch device has been implemented in place of an air brake 
from Ref. 9. 

The amount of torque that a brake can pass without slip- 
ping is a function of the type of brake and the net air pressure 
used to engage the brake. 

For a 35 inch diameter brake the torque is given by 

Q = 4550. p . (lb-ft) (25) 

x c *net 



where P net is the net air pressure activating the brake in 
psi. The brake glands which expand the brake shoes to make 
contact with the brake face expand inwards against the centrif- 
ugal force of the rotating axis. Consequently the air pressure, 
(P ) , required to counteract this force subtracts from the 

supply pressure (P ) . Thus the net air pressure to the brake 

s 

glands is given by 



P . = P -P 
net s c 



(26) 



The supply pressure typically increases upon inflation 
at the rate of 5 psi/sec and decreases upon deflation at the 
rate of 30 psi/sec. Thus, for a typical supply pressure of 
150 psi and with t in seconds, the deflating supply pressure 
is given by ' 
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P = 150 - 30. t 
s 



( 27 ) 



and the inflation supply pressure is 



P = 5.t 
s 



(28) 



The pressure required to counteract the centrifugal force 
is given by 



P = 5 + (5.9 10" 5 ).N 2 
c p 



(29) 



where 



N = 60 . n 
P P 



(30) 



therefore 



P = 5 + (3.54 10 3 ) .n 2 
c p 



(31) 



Combining the above equations yields the following expressions 
for the torque supplied by the brake: 

During engaging 



Q„ =4550 ( 5t-5-3 . 54 10 3 .n 2 ) 
X CL p 



(32) 



and during disengaging 

Q CL = 4550 (145-30t-3.54 10 -3 .n p 2 ) 



(33) 



Then, from Newton's law, the torque equation for the motor- 



propeller system is given by: 

dN 



2"i = ~(L -Q^ -Q, 



Qp -Q fr -Q cl 



(34) 
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Figures 20 and 21 give the results of using the above 
mechanical brake as a dynamic braking device. 

In particular. Figure 20 shows the variation of ship's 
speed versus time and Figure 21 shows the variation of propel- 
ler speed versus time. Notice the small change in the slope 
of the curve during deceleration. This occurs due to the fact 
that there is a delay in the application of the braking device 
and also that the torque produced by the moving water through 
the blades of the peopeller is still quite large. At t Z 10 
sec the change of slope is due to the application of astern 
torque on the propeller shaft. 

Since the use of a mechanical brake alone as a dynamic 
braking device will be associated with maintenance and heat 
dissipation problems, a combination of the above mechanical 
brake and the electric power dissipation methods examined 
earlier will be attempted. 

This could result in the smoothing of some of the problems 
usually encountered with the stand alone application of each 
system. Reduction of the size of the braking resistors or 
reduction of the heat dissipation on the mechanical brake are 
anticipated benefits. 

For this case Eq. (34) is modified by the addition of 
Eq. (20) resulting in 




( 35 ) 



where 



I tot I m ot +I shft +I p 
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For the case where the mechanical brake is used in com- 
bination with the linear variable resistor no significant 
advantage is gained as is shown in Figures 22, 23, and 24, and 
the reduction of the size of the braking resistor is insig- 
nificant since most of the dynamic braking is accomplished 
by the variable braking resistor. 

In particular. Figure 22 shows the variation of propeller 
speed versus time. It can be seen that for 0 < t < 2 sec the 
deceleration of the propeller is primarily due to the variable 
resistor since from Eq. (32) it can be seen that there exists 
a time delay in the application of the mechanical brake. This 
deceleration of the propeller is almost the same as that shown 
in Figure 14 where only a linear variable resistor is used 
for dynamic braking. 

Figure 24 shows the torque developed by the mechanical 
brake during dynamic braking. As can be seen it is not sig- 
nificant as to cause any change in the dynamics of the system. 
Figure 23 shows the variation of ship's speed versus time. 

On the other hand, the combination of the mechanical 
brake and the fixed braking resistor offers significnat ad- 
vantages as far as time requirements are concerned. Thus the 
time to stop the propeller reduces from 28 seconds, as in the 
case where a fixed braking resistor was used as the only means 
of dynamic braking, to 6.8 seconds with the addition of the 
mechanical brake. Similarly, the time needed to stop the ship 
decreases from t » 35 sec to approximately 25 seconds. 

Figure 27 shows the torque developed by the air brake. 

It can be noticed that for the initialization of braking 
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action approximately 1.2 seconds are required. Also, when 
disengaging should occur, it does not occur immediately but 
some delay is encountered. 
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VII. DYNAMIC BRAKING USING THE INERTIA OF 
MOVING PARTS AS LOAD 

Due to problems of space requirements and heat dissipation 
when electric or mechanical power dissipation devices are used, 
an attempt was made to find a way to minimize or eliminate 
these problems. 

Considering first the effect of ship's motion on the 
propulsion system, it is relevant to note the various mechan- 
isms present in the electric drive during astern maneuvers. 

The first stage in going CRASH-STOP is with positive 
voltage at the generator and positive armature current, the 
current being reduced at a fast rate. Since we want to apply 
astern conditions as fast as possible, the way to stop the 
propeller is to eventually reverse the armature current. 

The second stage begins the moment that the armature 
current reaches a negative value; that is, from this moment 
on the generator starts to motor from the propeller's kinetic 
energy 

During this second stage Eq. (10) becomes 

2n( VV Q E *Vn (3W 

where the relation for Q gen is given by Eq. (11) and for Q fr 
by Eq. (12) . 

As can be seen from the above equation the gas turbine 
will start to accelerate unless a way is found to reduce the 
effect of the reversed operation of the propulsion system. 
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Fig. lie? - Diagrammatic representation of the propulsion system 
dvring CRASH- STOP manouvres . (a) -First stage , 

(b) -second stage , (c) -third stace . 
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One way would be, if possible, to disconnect the gas 
turbine from the generator set so that overspeeding of the 
generator would not have any effect on the gas turbine. But 
in the model system we have selected, the gas turbine is 

l 

rigidly connected to the generator set, so this cannot be 
done . 

Because of the limitations on the overspeeding of the 
gas turbine given in Ref. 5 other ways have to be found. 

The limitations are: 

Power Turbine Speed 

Maximum continuous speed 3600 rpm 

Topping governor reference speed 3744 rpm 

Overspeed trip setting 3960 rpm 

Since maximum opposing torque is needed on the propeller 
shaft, this means that maximum current will flow in the 
armature circuit. This can only be accomplished by an effec- 
tive field generator current control according to the relation 

E = E - 33 (1) Volts (37) 

p- m 



But since 



E = K . i. (t) .N . 
g g fg gt 



(38) 



substituting Eq. (38) into Eq. (37) and solving for gives 

-r / _ E„ -33 

J fg (t) " i\r . .K (39) 



1 



This has been calculated in Appendix A. 
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Since the armature current (I ) will be maintained at its 

3 . 

maximum rated value, it can be seen from Eq. (36) that the 
only contribution that can reduce the gas turbine overspeeding 

effects is the term Q„, i.C., the gas turbine output torque. 

.b 

Thus, an attempt will be made during the second stage to 
reduce the gas turbine output torque. Since is a function 
of the fuel flow rate (W f ) and the gas turbine speed (N fc ) 
the reduction of the fuel flow rate in the gas turbine was 
selected as the means of doing this. 

At the CRASH-STOP command, the fuel flow rate is reduced 
to that corresponding to the idle speed of the gas turbine 
and the armature current is controlled and held at its max- 
imum value by the field generator current, according to Eq. 
(33) . 

Then, according to Newton's Law, the torque equation on 
the propeller shaft becomes 

- <W - Q p - Q fr <40) 

This permits rapid deceleration of the propeller. It must 
be noted here that in a real system the overspeeding effects 
of the second stage would be additionally lessened by the 
brush friction and generator windage and friction torques 
which at this stage have been assumed negligible. 

The third stage begins as soon as the generator voltage 
is reversed. From Eq . (38) this will happed when the propel- 

ler speed is at approximately 13 rpm; then astern motion of 
the shaft can be effectively applied. This is accomplished 
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by reversing the generator field current gradually, according 
to Eq. (39) , so that no overloading of the armature circuit 
will occur. 

Concurrently at this stage, the fuel flow rate at the 
gas turbine is increased gradually to cause the turbine speed 
to assume its normal operating value (3600 rpm) . 

When ship's speed approaches zero, a gradual decrease 
of the generator field current will bring the propeller shaft 
to rest. This can be accomplished either manually by the 
operator or automatically by the control system. 

The results of this method are shown in Figures 28, 29, 

30 and 31. In particular, Figure 28 shows the variation of 
propeller speed during the CRASH-STOP maneuver. The time 
to stop the propeller was found to be 1.5 seconds, i.e., in 
approximately 6 revolutions of the propeller shaft. In 
Figure 29 a plot of the law of variation of the field current 
is shown, according to Eq. (39) in order to have maximum 
armature current in the loop resulting in maximum deceleration 
conditions. Figure 30 shows the speed of the gas turbine 
during dynamic braking, where the overspeeding effects of 
this method can be noticed. 

The gas turbine speed increases up to 3660 rpm during 
2 

the second phase. Then phase three is initiated and astern 
conditions are applied. 



2 



In accordance with overspeeding limits given in Ref. 



5. 
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In Figure 30 it can be noticed that due to the fact that 
the fuel flow rate in the gas turbine cannot be increased 
immediately to the value corresponding to normal operating 
conditions, deceleration of the gas turbine speed takes place 
until the flow rate of fuel builds up to the required amount, 
where the gas turbine speed is brought back to normal operating 
conditions . 

Finally, Figure 31 shows the variation of ship velocity 
versus time. The time required to stop the ship, from a 
maximum speed of 30 knots, was found to be 25.2 seconds, 
which is quite compatible with the case examined earlier 
where a variable linear resistor was used as a braking device 
as is shown in Figure 13. 



( 
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Table II. Comparison of Results of Dynamic Braking Methods 
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VIII. DISCUSSION OF RESULTS 



The results to be discussed are shown in Table II. 

The use of a variable resistor as a braking device offers 
the advantage of fast deceleration of the propeller shaft in 
1.7 seconds from an initial speed of 230 rpm to 13 rpm. This 
amounts to stopping the propeller shaft in approximately six 
revolutions . 

Along with this advantage are disadvantages, however. 
First, the use of a variable resistor can only be associated 
with a combination of series or parallel resistors along with 
a voltage controlled timing mechanism. Large circuit breakers 
will also have to be incorporated. Space limitations and heat 
dissipation problems are additional disadvantages which will 
add to the complexity of the system. 

On the other hand, the use of a fixed resistor eases the 
problems of complexity and high maintenance costs associated 
with the timing mechanism and the circuit breakers, but still 
has the additional space limitations and heat dissipation 
problems though to a lesser degree due to the reduced size of 
the resistor bank (fixed) . But the decrease in complexity is 
paid in an increase of time results; that is the time to stop 
the propeller is increased to 27 seconds and the time to stop 
the ship from maximum speed is much greater than 35 seconds. 

The introduction of a mechanical brake gives almost the 
same timing as the variable resistor braking device does but 
still space requirements and heat dissipation problems are 



present. 
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As is shown in Table II, the combination of the variable 
braking resistor and the mechanical brake does not offer any 
significant advantage in time requirements and it almost 
doubles the space needed for accommodation of the extra mechan- 
ical brake. Similarly the reduction of the size of the resistor 
is insignificant. 

On the other hand, there are advantages to using a fixed 
resistor combined with the mechanical brake. There the time 
to stop the propeller is found to be 6.8 seconds instead of 
28 seconds when the fixed resistor is used alone or 9.9 sec- 
onds when the mechanical brake is used alone. The time re- 
quired to stop the ship from maximum ahead speed is found to 
be approximately 25 seconds instead of approximately 50 sec- 
onds when the fixed resistor is used alone, or 25 seconds 
when the air brake is used alone. But still, in this case, 
the disadvantages resulting from the use of heat dissipating 
elements still exist along with increased maintenance problems 
due to the addition of the mechanical brake. 

For the last case of dynamic braking examined, i.e., 
using the inertia of moving parts of the power turbine and 
the generator rotor, some assumptions were made. 

First, it was assumed that the fuel flow rate in the gas 
turbine could be reduced from the operating condition, to that 
corresponding to the idle speed of the gas turbine, instantan- 
eously. Second, complete controllability of the generator 
field current was assumed, so that constant maximum reverse 
armature current can be maintained. 
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Under the above assumptions, the simulation of the system 
gave the following results: 

Time required to stop the propeller from maximum speed 
to zero rpm equals 1.5 seconds. Time required to stop the 
ship from 30 knots to zero equals 23.2 seconds. 

For the operation of the system no extra equipment is 
needed such as braking resistors or brakes, nor any heat dis- 
sipating elements in excess of those already existing. The 
cost of the installation is not increased due to addition of 
extra braking devices and the cost of maintenance remains 
almost the same . 

In order to be able to compare the systems discussed, 
as far as performance is concerned, another assumption which 
was made for all cases examined has to be taken into account; 
that is, it was assumed that astern conditions could be 
applied immediately at maximum operating conditions. The 
effect of this assumption is that the time required to stop 
the ship will be increased by an equal amount for all cases 
examined. 

Based on the above discussion and having in mind that 
the scope of this thesis was to investigate the existing types 
of dynamic braking and specifically try to minimize and possibly 
eliminate the need for dynamic braking resistors, the last 
system was selected for further investigation. 

A block diagram of the open loop system in shown in 
Figure 5. 
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Fig. 5- Open loop system diagram 




IX. CONTROL SYSTEM DESCRIPTION 



A . GENERAL 

The system simulation was set up for the individual sub- 
systems outlined below: 

1. Gas-turbine-generator system 

2. Electric transmission system 

3. Generator field current control system 

4. Propeller-hull system 

5. DC motor. 

Subsystems 1, 2, 4, and 5 have already been discussed in the 
previous sections and there is no need to describe them here 
again. For references purposes the equations governing each 
system will be repeated here. 

Gas turbine-generator system 

The gas turbine-generator set speed is controlled by the 
influence of the engine torque and the generator load torque 
and friction, or 

M IPH -v Q fr 



where 



I =I_4- + I 

gt gen 



where the gas turbine drive torque is a function of the fuel 
flow rate and the speed, or 

q e = f ( "f ’ V 
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The generator torque is given by 



Q =K . I- . I 
gen g fg a 



and the friction torque is made up of brush friction and 
windage components assumed to be of the form 

Q : Q r K . N 2 
x fr br w g 



The engine is set to run at constant speed. 
Electrical transmission system 
The transmission equations are 



E 

S 






E 

m 





Ignoring the effects of armature inductance for maneuvering 

conditions, the armature current is given by 

x _ Eg - Em 
a R^ 



where 



R = R + R 
t g m 



Propeller and hull system and motor 

The propeller dynamics are described by the following 
equations : 



2n.! p 
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where 



and 



I = 
P 



I + I 
m shft 



+ I 

P 







I 

a 



The load torque on the propeller is defined by 



Q 



P 





.A 



where D is the propeller diameter, p is the local sea water 
density and 

A = (V 2 + N 2 .D 2 ) 
a p 



The coefficient C is defined by graphical data against 
the modified advance coefficient S where 



S = V / A' 

3 . 



1/2 



The ship hull dynamic effects are given by 




where 



and 



Tp = C,p . p .D 2 .A 
Qj, = f(S) 



and T^ is the ship's resistance as a function of the ship's 
velocity. 
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B. GENERATOR FIELD CONTROL SYSTEM 



The generator field coil equation is given by 



L 



'fg dt 



= V (t) -I 



c 



.R- 

fg fg 



where 



L,. is the field circuit inductance 

fg 



V (t) is the voltage of the field circuit 



R,. is the field circuit resistance. 

fg 



A closed loop system to control shaft speed using the 
generator field current to achieve speed changes would be 
appropriate for the following reasons: 

1. Speed control with a percentage accuracy better than 
the percentage system losses is not necessary and the operator 
could compensate for the losses if required. 

2. During turns using rudder control it is not desirable 
to compensate for the resulting shaft speed changes. 

3. The prime mover has been designed to run as a 
constant speed machine . 

The control system then is basically a field current 
controller. 

C. DEVELOPMENT OF THE FIELD CURRENT CONTROLLER 

From the discussion of the open loop system, it has been 
observed that the armature current must be constant in order 
to apply maximum opposing torque on the shaft. From the re- 
lation 



E = E - 33 Volts 
g m 



substituting the equation for E and solving for I_ gives 

g t g 
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I- = (E - 33) /K .N . 

f g m g gt 



Due to advances in solid state electronics and micro- 
processor chips the use of a microprocessor based controller 
was assumed. 

Comparing a microprocessor based controller with the up 
to now commonly encountered field current controllers incor- 
porating thyristor devices [7] , we obtain the following: 

1 . Control Room Space Considerations 

Because of the high density packaging used in micro- 
processor based controllers, large amounts of computational 
and control components can be placed in a relatively small 
enclosure to reduce control room space requirements. 

2 . Component Count 

Microprocessor controllers use "time* shared" com- 
ponents and thu S reduce component count. The same circuit 
that compute a 3-mode algorithm is also used to output out- 
of-limit alarms and extract square roots. One microprocessor 
controller, which can be used in lieu of up to six analog 
controllers and numerous computing relays, operates at least 
four final control elements. 

3 . Accuracy 

There is no question that digital devices are more 
accurate than analog devices. Using a 15-bit word length 
microprocessor uncertainty in any addition, multiplication, 
division or square root is about 0.003%. Digital devices do 
not suffer from drift problems. And with proper design, the 



57 



noise immunity of digital devices can be made superior to that 
of analog devices. 

4 . Reliability 

The reliability of a microprocessor based digital 
controller can be made much higher than the conventional control 
equipment that it replaces. High temperatures are avoided in 
microprocessor controller designs by virtue of the low power 
consumption of solid state electronics. Although the reliability 
of a single microprocessor controller may or may not be super- 
ior to a single analog controller, there is no question that 
a single microprocessor controller is more reliable than a 
system consisting of, let us say, six analog controllers, 30 
computing relays, and a panel full of logic circuitry. Use 
of a microprocessor also allows a number of "tricks" to be 
played by the programmer to enhance basic reliability even 
further. 

Self-test features can easily be incorporated into 
the program and failure of a self-test can be used to alert 
the operator that a replacement should occur. 

5 . Flexibility 

With analog and discrete components in a control 
system a change in control strategy normally requires a change 
in control equipment; these changes can be easily implemented 
by software changes in a microprocessor based controller. 

6. Cost 

Generally speaking for a single application the 
cost of analog devices may be smaller; but for a number of 
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applications the pertinent cost of the microprocessor based 
controller decreases rapidly. 

A simplified diagram of a microprocessor based controlled 
is shown in Figure 6 . 




Figure 6. A real time ^-controller system 
A dedicated microcontroller receives information from 
the input sensors and it is connected to the plant so that 
it may exercise direct control. 

Input to a microcontroller can be from: 
push buttons 
limit switches 
voltage sensors 
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temperature sensors 
pressure sensors, etc. 

The output of a microcontroller can be directed to: 
starters 
solenoid valves 
stepping motors 
displays 

servomotors, etc. 

Figure 7 shows a typical programmable controller processor 
diagram along with the associated units. 




Figure 7. The programmable controller/processor system 

The sensing devices of the process are connected through 
converters to the microcontroller. These converters (trans- 
ducers) reduce the size of the actual signal and then the 
signal is passed through A/D converters which change the 
analog signals into digital. After processing, the output 
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data are passed through a D/A converter which transforms them 
into analog signals to be used as control commands. 

In the actual system to be used the following inputs are 
required: 

a. gas turbine-generator speed 

b. fuel flow rate in gas turbine 

c. generator field current 

d. propeller speed 

e. field motor current. 

In the computer memory the following constants will be 
stored: 

a. generator constant K_ 

b. motor constant K 

m 

c. constant of voltage difference N = 33 

The required outputs depend on the control system to be used. 

For an effective control of the generator field current, 
it was decided that a positioning servo would be used to vary 
the field generator resistance. 



D. DESIGN CONSIDERATIONS OF THE POSITIONING CONTROL 
SERVOMECHANISM 

A positioning servomechanism can be designed as a 
second order system [13 and 14] , as shown in Figure 8 with 
an open loop transfer function 



G ( s) 



R(s) 

R(s) 



Wn 2 

s (s+2 jw n ) 



The requirements for such a system are: 



a. fast settling time 

b. restricted overshoot 
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c. almost critical damping. 

Selecting for our system 

settling time t = 0.01 sec 

maximum overshoot M =10% 

P 

and J = 0.9 we get from the relation 




solving for 

w = 444.5 rad/sec 
n 

Therefore the open loop transfer function becomes 

- R(s) _ 197527 

GU:> ‘ TTsT " s(s+&00) 

The closed loop system is shown in Figure 8 . 



1 *C8) ^ 

.yv i 


19759 ? 


c(s) 
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JTsTgoo) ' | 







Figure 8. Closed loop positioning servo 



The transfer function of the closed loop system will be 



R(s) . 

rTsT" 



w_ 

2 2 
s +2nw +w 
J n n 



19.76 10 4 
s 2 +800s+19.76 10 4 



E. FINAL REALIZATION OF THE jf -CONTROLLER ? 

The microprocessor based controller is shown in Figure 
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St1&orJ 




Figure 9. Condition monitoring and control of propulsion plant 

The output of the system will be 

a. output 1 - command to reduce fuel flow rate 

in gas turbine 

b. output 2 - value of calculated field generator 

resistance for maximum deceleration. 
Outputs are also provided for alarm purposes or self-test 
diagnostics . 

The operation of the system will be as follows: 

Upon issuing the command CRASH-STOP the controller 
is enabled for CRASH-STOP maneuvers. The necessary field 
generator current is calculated from Eq. (39) and then the 
required resistance value from the relation 

Rf = 

Rf If 

Simultaneously, the W p = idle is issued to a solenoid valve 
to cut the fuel flow rate in the gas turbine. Feedback is 
provided so that a false or slow response due to false operation 
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can be identified and corrected if possible. When the pro- 
peller speed reaches 13 rmp the system returns to normal astern 
maneuvers by giving the command to increase the fuel flow rate 
in the gas turbine. 

The above microprocessor based controller has been imple- 
mented for the CRASH-STOP maneuvers only but there is no reason 
why such a system could not be implemented for the complete 
control of the propulsion system. 

A complete diagram of the closed loop system is shown 
in Figure 10. The results of the closed loop simulation are 
shown in Figures 32 to 35. 

In particular, Figure 32 shows the variation of ship's 
speed versus time. Comparing this figure with Figure 31 
where the open loop system velocity variation is plotted, it 
can be seen that the time to achieve zero velocity with the 
closed loop system has been increased to 24.2 seconds, i.e., 
an increase of 4% more time than in the open loop system. 

This is to be expected since the response of the system takes 
into account the time response of the positioning servo- 
mechanism involved for the regulation of the generator field 
current. 

Figure 33 shows the variation of propeller speed versus 
time. A comparison with the results shown in Figure 28 for 
the open loop system shows a slight difference in the time 
required to reduce propeller speed from maximum to zero. To 
be specific, approximately a 2% increase in time is required 
for the closed loop system. 
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| Figure 34 shows the response of the positioning servo- 

1 mechanism where R c is the commanded position from the micro- 
1 controller and R^, is the actual value of the field circuit 
resistance . 

Finally Figure 35 gives the closed loop gas turbine speed 
response where it can be seen that the increase in speed is 
up to 3650 rpm, i.e., an increase of 1.5% above the rated 
speed in contrast with 2% increase in the open loop system. 
This is due to the fact that the small time delay of the 
positioning system causes a slight decrease in the armature 
current. 

An attempt to use an operator controlled open loop system 
for the CRASH-STOP maneuvers is considered highly inefficient 
since the operator's response to fast changes in system dy- 
namics will be very slow. On the other hand, for normal 
operating conditions and changes of propulsion plant dynamics 
an operator could quite effectively control the system when 

time response is not a serious factor. 

I 

Thus the proposed closed loop system using a microproces- 

If 

sor based controller is considered to be the most efficient 
way for controlling the propulsion plant dynamics under any 
operating condition, since this controller could, as stated 
earlier, be implemented for control of the system under any 
operating condition without much greater effort for program- 
ming it. 
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Fig 10- Closed loop diagram of the propulsion system 
with a microprocessor based controler . 




X. THE ROLE OF SUPERCONDUCTORS 



The relative flexibility and simplicity accruing to elec- 
trical power systems, in comparison to mechanical systems, has 
long been recognized. Specifically, electrical power and 
transmission systems provide flexibility of installation, are 
easily suited to automation and have a high degree of dependa- 
bility [15] . 

Despite these advantages, and others that could be obtained, 
the full benefit of electrical power systems has not been 
realized in marine propulsion applications. To be sure elec- 
tric drive systems have been used in ship propulsion, however, 
such use has been mostly limited to small or intermediate size 
propulsion plants . The primary reason that electric drive 
systems have not been adopted for large power plants has been 
that the weight and space for such systems is greater than 
that required for geared propulsion systems. 

As mentioned above, an electric propulsion system affords 
greater flexibility of installation. In contrast to a conven- 
tional geared propulsion system extensive shafting is not 
required and the location of the power source is not as re- 
stricted. In addition, the large noisy reduction gear could 
be eliminated. The noise generated by this component is of 
increasing military significance. 

For some time now it has been realized that when the 
temperature of an electrical conductor is reduced to a few 
degrees absolute, the resistance of the conductor decreases 
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to an immeasurable amount. The implications of this fact have 
considerable value for use in electrical machines. 

If the resistance in the windings of an electrical gener- 
ator or motor is removed, much greater current densities can 
be achieved with an accompanying increase in the magnetic flux 
densities within the machine. If this increase in flux density 
is obtained, greater power can be generated within a given 
volume. Thus, it follows that the volume of an electric 
machine can be reduced if superconducting windings are used. 

Superconducting machines also have disadvantages, how- 
ever. In addition to the requirement that a superconductor 
be maintained at a low temperature, there is a limitation on 
the amount of current carried that is exposed to a magnetic 
field. If the magnetic flux density increases then the per- 
missible current density is reduced. 

Similarly, the very fact that a superconducting machine 
develops extremely intense magnetic fields invites yet another 
complication. Since fields of the order of several kilogauss 
can be expected within the machine , disturbance of equipment 
in the vicinity of the machine will certainly result, as well 
as unbalanced loads on the superconducting field winding. 
Therefore it is necessary that a shield be provided to confine 
the magnetic fields within the machine. 

From Ref. 15 it was found that if a superconducting motor 
were to be used for the model ship derived in this study, 
certain variations on the values already derived in Appendix 
A would be necessary. 
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Table III gives a comparison of the values used on the 
model ship and those that could be used if superconducting 
motors were installed. The values for the superconducting 
machinery were extrapolated from Ref. 15, Table 3-1. These 
values are the results of a computer optimization method. 

For a generator set no data were available from Ref. 15 but 
it was mentioned that the diameter of the machine decreases 
rapidly as speed increases . For the model ship the generator 
rotor was selected to be 4.9 feet in diameter which is reason- 
able and compatible with that of a superconducting machine. 

This study on superconducting machines was made in order 
to be able to compare the values of the inertia used in the 
model study with another model. 

It is the author's belief that if superconducting machines 
were to be used instead of the conventional DC machines, then 
the variation of the specified parameters in Appendix A would 
be minimal and the results would not differ much from those 
obtained in this thesis. 

If, on the other hand, conventional DC machines were to 
be used, modification of the calculated parameters most probably 
would be necessary, since all calculations were performed on 
data extrapolated from different machines. 

It appears that the results obtained in this thesis for 
a conventional electric propulsion system, would remain sub- 
stantially correct if superconducting machines were used 
instead. 



69 



Table III. Comparison Between Conventional and Superconducting 
Machines 





SUPERCONDUCTING 


CONVENTIONAL 
MACHINE USED IN 
MODEL 


Operating voltage 


300 volts 


400 volts 


Operating current 


50,000 amps 


33,000 amps 


Motor speed 


200 rpm 


230 rpm 


Approximate diameter 


6.39 ft 


6.5 ft 


Weight 


113,000 lbs 


90,000 lbs 
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XI. CONCLUSION 



This thesis has presented a control model for dynamic 
braking of ships. 

The use of different braking devices, such as a resistor 
bank or a mechanical brake has been avoided by the use of the 
inertia of moving parts of the gas-turbine-generator set as 
load on the motor, which during the dynamic braking phase 
effectively operates as a generator. 

It has also been shown that by effective control of the 
generator field current, maximum deceleration of the propeller 
shaft can be achieved without exceeding the overspeeding limits 
of the gas turbine. 

The model used permits the study of different motors and 
generators as long as their characteristics are specified. 

It also permits the evaluation of a given propulsion plant 
with different types of propellers if their characteristics 
are also known. 



( 
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XII. RECOMMENDATIONS 



The control method that has been presented for the auto- 
matic control of ship deceleration during emergency conditions 
provides a realistic tool for further studies. The following 
are topics suitable for future investigation: 

1. Design of a microprocessor based controller for all 
phases of propulsion. 

2. Effect of sea-state conditions on the operation of 
the propulsion plant during CRASH-ASTERN maneuvers. 

3. Development of a complete propulsion system, in- 
cluding gas turbine dynamics and control based on microprocessor 
based controllers. 

4. Utilization of superconducting electric machines in 
propulsion systems . 
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Fig. 11-Valuos of the theoretical braking resistor 
for maximum deceleration. 
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Fig . 12-Approximated values of the braking resistor 
for maximum deceleration. 
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I 




:'xOt of ship^s velocity vs. time for the cases 
where the theoretical and the linearized resistor 
are used for dynamic braking. 
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Curve 1-Using the theoretical resistor' 
Curve 2-Using the line arized resistor 
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Fig. 14- Plot « of the propeller speed, vs. time during dynamic 
braking using the theoretical and the linearized 
braking resistors. 
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Fig. 15-Power dissipation diagram-Theoretical resistor 



PMAX=Maximum power dissipated 

PRB =Power dissipated in braking resistor RB 
PRM =Power dissipated in internal resistor RK 
PERC1= Percentage of dissipated power in RB 
PERC2=percentage of dissipated power in RM 
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T : ‘t. lR-Power dissipation diagram-Linear resistor 



PMAX=Maximum poser dissipated 
PRB =Power dissipated in braking resistor RB 
PRM =Power dissipated in internal resistance RM 
PERCl=Percentage of dissipated power in RB 
PERC2=percentage of dissipated power in RM 
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Fig. 18 -Plot of propeller speed vs. time when a fixed 
resistor is used for dynamic braking. 
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Fig. 19-Power dissipation diagram-Fixed resistor. 



PMAX=Maximum dissipated power 
PRB =Power dissipated in braking resistor RB 
PRM =Power dissipated in internal resistor RM 
PERCl=Percentage of dissipated power in RB 
PERC2=Percentage of dissipated power in RM 
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Fig. 20-Plot of ship's velocity vs. time when an air 
brake is used for dynamic braking. 
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Fig. 21-Plot of propeller speed vs. time when an air 
brake is used for dynamic braking. 
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Fig. 22-Plot of propeller speed vs. time when a combi- 
nation of a variable resistor and an air brake 
is used for dynamic braking. 



84 



Velocity 
(Knots; 



n 

L** 

i 



\ 



~T~ 

.in: 



\ 



> J 

Time, (.sec) 



n. 



'1 5 TRL ■ p 3 -• Du 
S'jT.RL" 6.02 



! ‘K 



u; > i 



'5 / INCH 
'5' INCH 



c iiN S3- L 
Plot h3- L 



Fig , 2 3 -Plot of ship ’ s velocity when a combination of 

frZ st n • ^ sl ftor and an air brake is used 
ror dynamic braking 

o 



85 






— t ■ 
• ■i 



Torque 

(ft-lb) 

xlO 5 



* ) 
l 1 

»n 



f 

/ 



A 




Time (sec) 


! h 


i r 


1 i 


J: 




i . it in. 2 : 


•^srnL c . -2.0Q 


U'OS'ItiCH 


Rl!N Nd- L 


vsniv r -c. 20 


l*: t S'Insh 


PlHT N5-3 



Fig. 24-Plot of 
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the torque developed by the air brake 
-Combination of the variable resistor. 
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25. - Fiot ot t*--. . • • / when an air brake 

and a fixed resistor are used in combination for 
dynamic braking 
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F»> 26 -Plot of ships velocity vs. time when a combination 
of an air brake and a fixed resistor is used for 
dynamic braking 
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Fi^. 2 f* a - Plot of the torque developed by the air brake vs.time 
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FI*. 2*. -Plot of propeller rpm s vs. time when the inertia of moving 
parts of the gas-turbine-generator set is used as load 
for dynamic braking 
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Fi*. 30- - Plot of g«e-turbine rpm s vs. time when the inertia of moving 
parts of the gas turbine generator set is used for load 
during dynamic braking 
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Fi* 31 . 



-Plot of ships velocity vs. time when the inertia of moving 
parts of the gas turbine-generator set is used as lodd 
during dynamic braking 
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-Closed Loop-Plot of ships velocity vs. time during dynamic 
braking and astern motion 
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F»^. 35 . -.Closed Loop-Plot of propeller rpm s vs. time during 
dynamic braking and astern motion of propulsion 
system 
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Fig . 34-Closed loop response-Variation of generator 
field resistance 
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Fig. 36-Closed loop response-Variation of the generator 
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APPENDIX A 



DETERMINATION OF ELECTRICAL AND MECHANICAL 
CHARACTERISTICS OF THE SYSTEM 

A. DETERMINATION OF THE ELCTRICAL CHARACTERISTICS OF 
ELECTRIC MOTOR AND GENERATOR 

Since data of an actual working system were not available 

at the time this study was conducted, an effort is made here 

to establish the electrical characteristics of the motor and 

generator, based on general known characteristics [4] . 

The maximum torque on the propeller due to the water 

moving through its blades can be found from the tabulated 

data of the propeller and the equation 

Q = C .p. D 3 (V 2 + (N.D) 2 ) (1? 

p q p 

Thus for the maximum number of shaft rpm's this was found to 
be equal to 

Q =583580 ft-lb (2) 

P 

Assuming the frictional torque on the shaft to be constant 
with a mean value of 22,500 ft lbs, then the maximum torque 
exerted by the motor according to the equation 

2y tl |p =Qm - Qp -Qfr (3) 

under steady conditions is found to be 

Qm = Qp + Qfr = 608580 ft-lb (4) 
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Translating this to horsepower gives 



105 - -Mo§- = 26650 • 87 H P 



( 5 ) 



Expressing this value in kilowatts we get 

P = Hp x .746 = 19881.55 Kw 



( 6 ) 



For marine applications the maximum operating voltage of 
motors is 600 V. Therefore the circulating maximum armature 
current will be 

T P( watts) „„r n /"7\ 

X a ~t v5It Z T 33135 - 9 “"P 217 ® 3 (7) 

By now motor efficiency of 1 has been assumed. 

For this size motor almost 6% of the power is dissipated 
in the armature resistance, so 

Rm + Rg = R t = = 0.00109 ohms (8) 



Assuming equal internal resistance for generator and motor 
we get 



Rm = Rg = 5 10 



-4 



ohms 



(9) 



The back emf of the generator becomes 

Eg = 600 + la . Rg = 616.5 Volts (10) 

and of the motor 

Em = 600 - la . Rm = 583.5 Volts (10a) 

under steady state conditions. 

So for ahead maneuvers at maximum speed 
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Eg = Em + 33 Volts 



( 11 ) 



and for astern maneuvers 

Eg = Em - 33 Volts 



(11a) 



The generator back emf on a D.C. machine is given by 

E = k .ip . n 



( 12 ) 



from which solving for K 

K = E /(i f .n) 



(13) 



Assuming, at steady state, a field current of 10% - 12% of 
the armature current 



then 



I_ =40 amperes 
f max ^ 



Therefore at maximum operating conditions at steady state 

Em = 583 volts 

N , = 60 rp/sec 
gt 

N = 230 rpm 
prop 

Eg = 616 volts 
So 

Kin = 0.06326 and Kg = 0.25666 



Similarly, from the torque equation of a D.C. machine 
Q = K .i f -i a 



(14) 



having in mind that the output torque of the motor must equal 
the sum of friction and propeller torques, is found to be 

0 = 608580 ft-lb 

miax 
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Solving Eq. (14) gives 



Km = 0.459153 

For the generator, considering maximum operating conditions 
and steady state at 3600 rpm the tabulated data give 

Qg = 38000 ft-lb 

therefore from Eq. (14) 

•Kg’ = 0.0286697 



The above results are summarized in Table IV. 

B. CALCULATION OF PROPELLER MOMENT OF INERTIA 
From [16], propeller weight is given by 

W =K.D 3 (MWR)(BTF) 

where D is the diameter of propeller in inches 

D =13.94 ft =167.28 " 

I 

I I 

MWR is mean width ratio = .8 for a three bladed propeller 

I ! 1 

i BRI is blade thickness fraction and is assumed to be .05 
| K is assumed to be .26 for three bladed propellers. 
Substitution yields 

W = 48681.66 lb 
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Table IV. Summary of Characteristics of Motor and Generator 



AA 


DESCRIPTION 


MOTOR 


GENERATOR 


1 


Torque constant 


Km'=0. 459153 


KG*=0. 0286697 


2 


Electrical constant 


Km=0 . 06326 


KG=0. 0042916 


3 


Back emf (max) 


Em=583.5 volts 
-4 


Eg=616.5 volts 
-4 


4 


Internal resistance 


Rm=5 10 ohms 


Rg=5 10 ohms 


5 


Internal inductance 


negligible 


negligible 


6 


Field voltage 


E_ =400 volts 
fm 


E£ =400 volts 

fg 


7 


Field resistance 


R ^=10 ohms min 
mf 


R_ =10 ohms min 

fg 


8 


Field inductance 


L £=20 hrs 
mf 


Ij£ =20 hrs 

fg 



At steady maximum field current I f = 40 amps 

state 

maximum armature current I = 33136 amps 

Notice that the difference in Eg and Em is 33 volts; therefore, 
without overloading the armature circuit, application of astern 
motion starts when Em = 33 volts. 
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Therefore propeller mass 

m=Y =1511 ' 85 lb-sec 2 /ft 

The radius of gyration is assumed to be .25 D so 

I ' = ( . 25D) 2 m = 18361.7 ft-lb sec 2 

An allowance of 25% increase is allowed for the inertia of 

entrained water in the propeller. Therefore 

I = I (1+0.25) = 22952.2 ft-lb sec 2 
P P 

C. COMPUTATION OF MOMENT OF INERTIA OF MOTOR ROTOR 

Since no data were available for a D.C. motor this size, 

considering data which apply for a synchronous motor of 15000 

Hp rating output power and which give the weight of the 

rotor to be approximately 69,000 lbs, a rotor weight of 90,000 

lbs was selected [12] . Therefore 

m = — =2795.03 lb sec 2 /ft 
g 



and 

I =1/2 m r 0 2 
m i 

For an outside diameter of the rotor of ( 2 m ) 6 . 5 ft we get 

I = 15042.72 ft-lb sec 2 
m 

D. COMPUTATION OF MOMENT OF INERTIA OF PROPELLER SHAFT 

Assuming a shaft length of 50 ft with outside diameter 
D2 = 1.6 ft and inside diameter D1 = 1.37 ft gives: [12] 

Volume = 2n(r2~R^)L =36.12 ft 2 
104 






The specific weight of steel is 490 lb/ft 

w = 17702 lb 
m = 549 lb 



so weight of a shaft 
so mass 
and 

Total moment of inertia of motor-propeller system 



I = 175.9 ft-lb sec 
s 



I = Ip + Im + Is = 38171 ft-lb sec 



E. CALCULATION OF GENERATOR ROTOR MOMENT OF INERTIA 

As in the case of the motor and taking into account the 
maximum number of revolutions of the generator a rotor weight 
of 40,000 lbs [12] was assumed 
Therefore 

m=— =1242.24 lb 
S 

and assuming an outside diameter of D = 4.9 ft we get 
Ig = 1/2 m R 2 = 3760.7 ft-lb sec 2 



( 



F. CALCULATION OF TURBINE ROTOR INERTIA 

From [9] the following data were extracted 
Low pressure compressor Jp = 586 

High pressure compressor Jp = 489 

Free turbine rotor Jp = 5009 

where Jp = polar moment of inertia. Therefore 

1 =-§-.2 'IB = 155 - 56 ■ se ° 2 
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Therefore total moment of inertia of gas turbine-generator 



system 

Igt = Ig + Igen = 3916.243 ft- lb sec^ 



G. CALCULATION OF THE WEIGHT OF BRAKING RESISTORS 

In order to be able to select the system that will be 
applied some criteria about the size of the braking resistors 
are needed for each case [17] . 

The energy that has to be dissipated in the resistor is 
given by 

E = i a (t) in joules/sec 



The energy .which has to be absorbed is given by the area under 
the curve of the previous equation. 

If the resistor were of copper with specific heat of 
0.0918 BTU/lb°F, for an acceptable temperature rise of 100 °C 
the weight of the resistor will be given by the formula 

Weight in Ihs = Energy absorbed in joules 

Weight in Ihs o>0gl8 _^ x212 , p x 1064 8 joules 

The energy absorbed in joules by the braking resistor in 
each case was found by using Simson's discrete integration 
formula . 
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Table V. Weights of Necessary Braking Resistors 



Method of 


Energy Absorbed 


Weight of 


Dynamic Braking 


in Joules 


Resistors in lbs 


Use of a Variable 
Resistor 


1.38717 10 7 


675.7 


Use of a Fixed 
Resistor 


1.7621 10 6 


85.883 


Use of a Mechanical 
Brake with a Variable 
Resistor 


1.3870 10 7 


673.3 


Use of a Mechanical 
Brake with a Fixed 
Resistor 


1.2549 10 6 


61.131 
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V V >/ r 



’ V, 



r; tr.cn cost 
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*3. -Map of propeller torque characteristics 




Fig. 39-Propeller torque characteristics ,Cq 
as a function of the second modified 
advance coef f icient , s 
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4.- Map of propeller thrust characteristics 




Fig. 40-Locus of thrust coefficient Ct and second mo- 
dified advance coefficient. 



Ill 



5.- Map of ship’s r e s i stance , propuls ive 
coefficient and propeller speed vs. 
ship's speed. 




.72 

.70 

.68 

.66 

.64 

. 6 ? 

.60 



Fig 41-Plo. resistance (R) , propeller 

speed (Np) and propulsive coefficient 
(n ) versus ship's speed (V). 
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APPENDIX B - COMPUTER PROGRAMS USED 



v Tf-JS PRCCPAM WAS WRITTEN IN DSL/360 SIMULATION 

* LANGUAGE . TRANSLATION OF THE PROGRAM WOULD St 

k i ■ r - 1 H J I P E C TP ANOT HER SINUl ATiCA ■ LWN'"i. i Ar ,H — t; 

TITLE FPDGRAM A 



* 

* OPEN LCCP SYSTEM 

$ 

b p -l- t -gp m i n at i un cr- T . nr v a l u - e & -- g f the- t h or r i .. t g a i 

* variable resistor 



v " ••• 1 ■ 

$ 

$ 

=1 VARIOUS CASES OR DYNAMIC BRAKING ARE INCLUDED JN TH'S 

* PROGRAM , AND ARE THE FOLLOWING : 

■A 



* 

* 1 . -DYNAMIC BRAKING LSI MG A LINEAR VAR^A^LE 

| v RFiVSTCK . 

* 2.-DYNAiv*c BRAKING USING A F U ED RES 1ST NR 

_* 3 .-QYN'A-«-(. TEARING UjJjSiG A r 1 E C I i A N ■ C A L 

~ A’N TYPE bRAKE . 

* 4. -DYNAMIC BRAKING USING A COMP! N ATT ON' 

—a U£ — A VAR TABLE I.' F 3 I S TUP — AA1 Cl— T-l i-G . - M £ C H A N V . C, N I 

•> - A ( P BRAKE 

* 5. -DYNAMIC BRAKING USING A COMB! NATION 

tt rr^ mi A f< ti K A K t AND 3 rim*) !’ E ) 1 S T 0 w ; 



* WHAT THE USER HAS TO DO IS TO StLECT A METHOD BY 

* INDICATING IT BY T TS NUMBER ABOVE 

i P2j2_ P XA-M RL'E CASF=2 MEANS THAT OYMM'C — BR AK. T N C 

* WILL HE ACEOMPL'SHED BY THE USE OF A FTXEP 

* PES T STOP 






IN ADDITION 



* TM<= USER HAS TO SPEC T n Y T^E FOLLOWING : 

& ^A, ) ..---M.AX ' MIIM PERN X .S.S4 B-L-fe— C UR K FN T ,'N THE. ARMATURE 

* CIRCUIT (THIS IS ONLY USED IN' T Ht CAI.CU- 

* L AT I ON OF THE THEORETICAL VALUES OF THE. VARA APLE 

RES f ST OR- ) 

* IB) - VALUE? OF THE ELECTRICAL CONSTANTS CF THE 

* MCTCR AND GENERATOR, M AND KG RESPECTIVELY. 

1 r'C )" - VALUES l' L THE MECHANICAL CDNSTaNTS”OF' MOTOR 

* AMD GENERATOR, KMM AN C KGG RESPECTIVELY . 

i UIU— =- -VALUES OF. THE INTERNAL R. ES " SI/UCCE-S- ■ CF . MO LOR 

* AND GrNERATTR, RM AND RG 

* (E) - VALUES FOP. INERTIA CF PROP ELL FR , 

b SHAFT A NO MOT)K SYSTTM AMD WrttH^S P T: H — Tt+= 

* INERTIA CT GENERA TCR RCTCR. ,PPWfcP TURBINE SYSTEM 

* AS 'N APPENDIX A CF THE THESIS, 

lit “ 



* 

Tt ! T I. (> TRAP? 

intger nflot 

Intgek case 

CONST |\PIU! = 1 
I* 

A: g£ tcqm f VATTMiM OF HUPP'I 'FC PARAMETERS 

J. 

CLNST I-?ai7l. ,C=13.94 
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WTEPM, nation of shtp pa paveteps 



;CI\ST P=2. 73970+05 



PFTEP M I NAT ION GF MOTOR CH AH AC T L-R I ST I CS 
■ C 1*5 T ■ + 9 - 0 T ) 6 3 ? L r 1 K “P « Q . 4 b 9 1 f ? 3 , < . ; 1 - Q , Q Q Q £ 



PET ER V 7 n a t !QN O c GENERATOR CHARACTER? ST K S 



CNST KG- 0.7 56666 » K GO— 0 . Q 2 B 6 6 9 7 iRG = C.OOC5 



jl'-RTTA CF GAS-TUPBTME/GENEPATCH SET 
I f ' l - S - T VC- - l r ‘1.6 



OTHER CONSTANTS 



Ci\ S T P' f = 2. 14159, P = ?. 

mil LAi£=i 



in t t:ai. ccnottjcns 



NCCN £ C 1 = 30 . t T C 2 = 3 . 8 3 3 3 t I C3= 0 . 0 » I C 9 = 6 C 



ERIVAT7VF 



l = :mgp». ( LCJ .unrT i 



X= r MGKL ( T C 2 » U J 
NGT = 7. NT GRL ( IC4 , GTPOT ) 
•mw (v) 



Vs =T AFW( V ) 

V P= ( 1. . - lv ) * V 



T - T AN. * ( V > 

A = VP^*7 + (N-»C)**2 

B = S 0 P T (~ A ) t 



S=K':M!/B 
CT=TABCT{ S , VP ) 



TP= ( l.-T)*TA 
TPR = ( TP-P.T ) /m 



~ TJ P C' T = T p R 

V=U 

C G = T A B C Q ( s .VP) 

cp = c;g-p* (0**3 )*A 
N = ?. NTGOi. ( IC7 ,NDOT) 

" ' - V -V f P 



-m 



-E-N! * N p F A I.. 



E 0= X 0* I F G l * N G T 
(A=72256. 



rf e i.~ ( r ,, ~: : . > / ( kg -' -7 ' 0 t i 

QWG=- r A** C P*KVM 
PRP= ( V A ** 2 ) *P B 



P p N' - ( T A « * ? ) ^ R M 

P P A X = c M * I A 

_Uf) T — !- M-T Atov 



PfRCl=( PR8/PPAX )*\00. 
PFPC2= ( PRv/FMAX) 1 0 ") . 
tty -TV ‘ ■•E 



M\AP 1C 



~ 1 T ( *i . L T « 0 « 2 l 6 u ) - GG TO 700 
T F ( C.A S F . EG . 1 ) GG TO 100 
_I P ( CASF.EQ.2 ) GO TO 200 



THcASE.fcC .2 ) 
i F ( CA S F . EQ .9 ) 

_L F ( r a £ c ~ c r . 5 i 



GU TD BOO 
GO TO 300 
-GO TO 330 
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uc 

Fr& 



k 8= ( E M — I A* P F )/! A 
If l±l -I T.r-rnis ) 
GC TP 120 
R 8=0 . 0 0 1 5 

CC.F'T - NLJF 

GC. TO 600 



r.n to i i n 



un-M^r.T •«&TM G F N F R A T R AMQ AND A P.CR5TANT £ f c , 7 S T C P 
2C0 RP=0 . 0 l 7 
:fm=40. 

; t k C 

GVCT=0.0 

GC TC 6 OC 



St CF AR AIR CLUTCH 0 R A K 7 >N| G PEV.tCE 

PC-5 . +0 . ?! OAMN**?) 

P 5=5.* TIM 
QCL=4550.*( FS-PC) 

crtrr^rr o 

IF(QCL.LT.C.C) GO TO 30 L 

GC TO 302 

:0 L 001=0.0 
302 CONTINUE 

PNET-l 50 . 

0CLF=455C.*FRET 
JF(CCL.GT.CCLM) GC- TO 303 

PC T O "SOT 

3C3 CCL = OCU R- 

S04 POR T 7 MU F 

? F(CAGE.FC .A) GO TO 100 
1F(CASE.EC.5) GC TO 200 • 

GC t^ - 600 



t'0'0 SO^P=;juP T - C P-OF P.-OC l +0 MG +u M G'(' 

GO TO 900 

700 T E ( C A S F . G E . 0 . AN C . C A S F . 1. T . 6 ) CC-T0-7-UX 
GC1.=0 . 0 
C-C TP 71.1 

ttT j PC- 1 :, J . 2 1 2 A * ( N « * Z: ) 

K = T 7 ME— 6.7 
PSZ =PS-3C. * K 

lvUl. = 4 5 60 . * '( FRI-PI. ) 

iF(CCL.LT.C.O) GO TO 701 

or TP ,702 

7C1 GCL = o. 0 
702 CONTINUE 

T t t C’O'-O.O 

PP8=0. 0 
FRv=0.0 

n-RU=o'.'c — 

p EFC2 = 0 . 0 

U*.L . 0 

E F = 0 . 0 
RP. = C.O 

— 

I AST=— ?>3 P 0C . 

QAST=KMf'*"FM=TASr 

GT i f J 8 00 — 

8CC SIJMO=QAST-CP-CFR+CCl 

M |GLLLXAUJ£ 



■ 



CFF = f-f.MS!»(N ,-22 500., 0.0, *2 25 JG. J 

Ar nT=SU?«C/(?.*Pt* T ) 

(\k ; mi. =b 1 > • * r> 

rrki=pr 8/ icoooo. 

0 g M i - p S v / i r r r n . 

PVAXI = PN‘AX/10000 3. 

IAl=7A*C.C1 515 

nr, /i oo oo, 

CP2=GP / l C C C C . 

CCL l=OCL / ICOOCO . 



L£ 

ftl F I N H' M = ?. 5 . , C EL T = . I , D E L 5 = . 1 

T . I , FM , P R B , Vn. T , CP , V , T A , P RN t p E *c 1 , CP R , NF E A L , R B , PM AX , DP P C 2 , C M 0 , . . . 

I r r A-S r P '-A S - I QM - 3 - Ty , S -t ) r 5 , E G i C Q , P - S , G C - L . '■ ) 3. , , >T - |- , c T , p C . Q C, F A KG T . , . . . 

i F G 1 » W F , I.A,CE,VP,Vv,T,4,B,TPK,K,PSI 



THE FOL LOWING STATEHFNTS ARE APPLTCA81 P in THP 
PI T t t*\t. paCKAOH FHP pst S * U I j T , r N I A^nil^GF 

AT Ni.P.S . THE STATEMENTS SHOULD ti- CONVERTED IF 
CTHFP THAN THIS PLOTT T \ G PACKAGES ARP AVA'l_ABt c . 



CALL n«wC( l, l, TIME, V) 

C A i_ L DRUG ( 2 ,1 ,T ! me,NREAI ) 

mrr — i rwi 

CALL 3RW G ( 3 * 2 » T I N F , P RM .1 ) 



CALI. 1RWHI 


L 2-j-3.,T J mf , PERC l j 


1 


CALL ORWGI 


! "> ,4 ,T I HE, P^PXR 1 


\ 



CALC DPHG( ? ,5, TIME. P M A X l ) 

■ CALI — l)PWQ( -4 , 1. , TI ME , I AT i 

CALL ORViO ( V,? ,T \ ME, EM) 
CALL 0RWG(5,1,TIME,0CH) 

trat 

CALL E NO 3W ( \ PLOT ) 






TABULATION OF SPEEC REDUCTION COEFFICIENT VS.S°FED 



FUNCTION TAHU(V) 

t ) t ' * E fc S 1 0 M VT (13) , ' .IT( 1 .3) 

DATA VT/O.o, 2. 5, 5. ,7. 5, 10. .12. 5, 15. ,17. 5. 20. ,22. 5, 25. ,27. 5. 20./ 

DATA NT/4*C.C,.J05,«01,.32,.C-;3,.C45,.Q45,.008,.02,.OC4/ 

iHv.ij 1 .’30.) go irn 
!3FLV = 2. 5 

AaTFTXC V/nFI. VH-1 

SLCPW=( WT ( IS + 1 ) - W T ( N ) ) / ( V T ( N + 1 ) - V T ( N ) ) 

TABW=SLOPW* (V-VT ( M ) ) +WT(N) 

ir ETi JR *i 

T AHw=0 . 0 
RETURN 

tnh 



TABULATION OF SHIP’S RESISTANCE VS. SPEED 



__ FUACTTOM TARR(V) 

GTV-f SION VT (13 ) , k T T l ITT 



OAT A VT/O.C.2.5, 5., 7. 5, 10. ,12.5,15 
UlLA P TT / 2*0 .0, 7000 .. 13000 . .PEPPO - 



1120 TO 3 • , l 7 2 C 0 0 . ,225030. , 23 3000. / 
<F(p.GT. 3C.) GO TO l 



. ,17. 5, 20., 22.5,2 5. ,27.5,30./ 
,3 50 0 0 . ,5 70 00. ,BO O O C . , 1 C TCOC. 
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0FLV=2. 5 

K='i cT y (v/cf'.vi + i 

' b U 1 • 1 r - t T. TT ( | k ^ 1 1 - T T ( ' -11 1 / f V T 1 ,"1 + 1 ) - V T ( ‘ 'i ) ) 
TABP-SLHPP* (V-VT(M) ) +PTTIN ) 

rttioo 

TARP=2800CC. 

RETURN 





T A BUL AT 7. CM OF THRUST RFDUC T ? ON COEFFICIENT VS.SPEEC 



FCNC T !ON TART(V) 

— DT N FNS • ON V T (13) , TTNIO ) 

C A T A VT/O. 0,2. 5, 5. .7. 5, 10. ,12. 5, 15. ,17. 5. 20. ,22. 5,25. ,27. 5, 30./ 

CATA TT/5*C.Q».0l» .05, .07, .CP, .075, .072 ,. 16, .02/ 

gt:t.y- ? rg 

N =T F T x ( V/CELV ) +1 
iF( V . C. T .30 . ) on TCI l 

SlUPf = ( T T ( A + 1) - TT ( N: ) ) / { 7 T ( N + 1. )-V T ( U ) ) 

TABT = SI_OPT* (V-V T <N) ) +TT (N ) 

p Exuajj — — — — 

TAtiT-=0.0 

RETURN 

— tr-fl 



T A HU LA T T UR H F rTTRTTSl CPfcFr iC IEM VS.SLCPM iv[’.i fFTFC 
AO V A NIC r CO E FF If 1 ENT 



FUNCT* ON TAECT( S , VP ) 

D.TNFUS: ON ST (21 ) , C T ], ( ?U ,C T 2 



-ti- 



0 A T A ST/-l.,-.9*-.b,-.7»-.6»-.5»-.4»-._>* .7, 

1. A, .5,. 6,. 7 ..8, .0,1./ 

l! * T A C T ?/- . 4 , -. i 5 , - . 35 , . .>5 , . 1 , . 1 5 , . 2 , . 2o , . 36 , 



1.45,. 42, .35, .37 , .22, .33, .50/ 

-CAT.\ CTl/- .4,-.2. 2B,~. 35,- 



■8.. 



1-. 35, -.3 4, -.2 5,- .2,- .13 ,-. 05 , .Cl , . 10, .21, .45/ 
CEL c j= 0 • l 

■■ hi - T F ■ *< ( ( & t l. ) - /oc - LS - y r\ 

•F( VP .LT .0 .0 ) GO TO 2 

SLOCT 1 = ( C T 1 (N + U-CTl ( N ) 1 / l S T (N+ 1 ) -ST ( N ) ) 

I A R ( , I = bLLC ! 1 * ( S - S 1 (Mil + ( I Id' 1 ] 

GC TP 6 

-SCR. T?.= <CT 2 (N + l ) -CT2 UU 1 / ( ST I N+l > -ST( M ) ] 



TA8CT=SLQCT?.*( S-STIM ) ) + C T 2 ( M ) 
RETURN 

-E+rf) 



1,0. c, 

44Y.’3 0 

■^.■■36,- 



1 , . 2 , . 3 , 
. Z 1 , . 4 , 

K.-,7 c , 



I A BUL A I TlJN FTF rUTCTTF CUFFFTTTTETri VS.SHC'M. ^TTTTFTEX 
AOVANCF COEFFICIENT 



FUNCTION TAECCt S,VP) 

~ E T MR ' .ST ON G - T 'f- 2 ' 1 ) ,C01 ( 31 ) , CO 2 ( 21 ) 

CAT A ST/-~1., — .9, — .F' , — .7, — .6, — .5, — .4,— — .2, — .l,0.,.l,.^,.- , 

1.4, .5, . 0 , . 7 ,.4, .0,1 ./ 

Urn Lo.l / - . 04 , - . O^+p , - . up , - . , - . t e , - .uc-i , - , - . 'Jt*+ , - . c 5 1 ,- 

i-.037,-.055»-.05, - .035,-.U26,-.0l5,-.00d, .0C5 , .02, .002 , .C7/ 
■■DATA r 0 2 Lz -PR,- - 0 25 01 5 . . 0 1 . . 0? . . 22Jt+ . J4 ? . . Th . . 0 ~ . . C 4 5 . 

l.Cer>,.07,.C62,.C6,.C50,.056,.C5,.Cfc/ 

DELS=0 . 1 

-t T=ir ;x( ( 5 H. ) /DLLS ) n 

lP( VP.LT.O.C) GO TO 2 

SLCCOl = ( C G 1 ( N + l ) -COT (N) )/ 1ST ( N + 1 ) - S T ( N ) ) 

T AM.U=SCITC01* ( S- S I l N ) ) TUTJT (N 1 
GC T fJ 6 

—Si rro?= (ro?(N + i )-rc? (Ml i /( ST(N-t-i )- ST.LkLJ 

TABCO=SLOCC2 5 ! ! (S-ST (N) ) +CQ2 ( N ) 

RETURN 



TTT7 

CL5-* 
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tr.KD 



TABULATION DF PARAMETERS OF THE PHIN'C MCVE: 



FLM'T 7 ON T M'QfcM V3 ,V*F ) 

UTWEMS T GIJ VST (5 ) , V-;FT ( 10 ) , 0 E T ( IOC, ICO) 
V jT ( I) tQ 



V 3 T ( 2 ) = l G C C . 
V3T(3) =2 OC C . 

V ’ HR)- TWZ CT 
VST ( 5 ) =4000. 
UFTM ) = ? ? 0 C . 



R FT ( 2 ) = A30C , 
UFT ( 3) =53CC . 
UFT C - " 



hFT(5J»73CC. 
wFT (6) =3300, 



xf t rrr = t?n: 

WFT (3 
In FT ( 0 



) =10300. 

) = n sen. 



W F T ( l 
OFT (). 
LX 



J )=122CC. 

1 )=2C(!00. 



GET ( l , 24 ~2 FS 0,0 . 



CET( 1 
G 6 T ( 1. 
O c T - Cl - 
OFT (1. 
GET ( l 



3 ) = 
A ) = 

e> ) = 

7) = 



37CC0. 

A5CC0. 

ftp COO. 
7CCOO. 



~3T= 

,9) = 



7 7 C 0 0 « 
>3 ft 500. 



OFT ( 1 
QFT( l 



•T ( (. , 4.C) 000 , 



GET ( 2. 
GET (2, 
O c :TC?, 
OFT (T.i 
OFT f 



l) 
2 ) 



T=- 



125CC, 

:?C000. 



A) 

ft ) 



ill;. 
2 A 500 . 
■ 7 T 5 0 0 . 



OFT ( 2, 
OFT ( 2 , 
■ Q . ET I ? .. 



6) 
7 ) 



= AECCO, 
: 5400C. 



) ^- 5 v OO-: 



OET (2, 

OFT (2 



OFT ( a. 
Q FT (ft 
QFT (3 



9 ) = 

10 ) 
tr^ 
2 ) = 
3) = 



64000. 

= ftFCCO , 
3 3 C C . — 
13C00. 
It 500, 



OFT ( : 
OFT ( p 



2 A 0 C 0 , 
29000. 



4T^ 

5 ) - 



■ GET(';,ft).= 3ACCa. 



OET ( 3. 
GET (3, 



7) = 
3) = 



? C 5CC. 

45CGO. 



t jftT( 3, - 7) *QCCCO. 



OFT (2 
WET (A 



10 ) 
l ) = 



-55000, 

5CCC. 



‘Lb HA, 

G'FT (A, 



TT= 

3 ) = 



9 C C 0 . 
L2C00. 



-C E T ( 4 , 4 ) = 1 7 c -:o , 



GET (4 
GE T ( A 
-g' r T(-V 
QF T ( A 
GET ( A 



5 ) = 
ft ) = 



) -2 2 ' 
3) = 

9) = 



2 2000 . 

26000. 



3 ACCO. 
32500. 



l^TT 
GET (ft 



TTTS 

1 ) = 



=42000. 

3 CCO . 



- 4 r T(ft , 2) =6 *00 , 



GET ( ft, 

C- 



A ) = 



•fr F'T ( ft , 5 ) - 1 
6 ) = 

. 7 ) = 



ICCCO, 

12000. 



GET (ft 
OET (ft 



UT-' T ( 5 , 
GFT (5, 
-ft ft i 5 



TTT= 
9 ) = 
-UU 



2CC00, 
2 A COO. 



t mrcrr 

3C5C0 . 
s 24000 . 



iF( V3 
i F ( V 3 , 



T 

GT 



.0.0 GO TG l 
.ACCO.) GC TO 2 
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X 3 -\'3 

gc: 3 n too 

I Ai = ' . ) 7 

or th ioo 

X 3 = 400 0 . 

GO rn ioo 

NO IF(mF.LT.32GO.) GO tp 3 

730 3 . ) GG - TO / , 

v» F X — w p 
GC TU 0.00 

AFt-33 TO. 

cc rn 300 

■*FL = 12 300. : 

200 1 = 1 FIX ( X 3 / 1 0 0 0 . ) +1 

J = T F ’ X ( (WFl-?30Q.)/1000. )+l 

OJ g- y - G-v - 3 - T . ( t ) — 

L V- = V- FI-V./FT { J ) 

UFLOR= <np/ 1000. )*(OCT( T +1 , J )-QET ( 7 * J ) > 

r o- r i v -c t m/tcco. rrrnrm . .if b rt -v, j) > 

TAHOr=OET ( T , J ) +BELOR +DELGh 

R c ThRN 

END 



THE FQLLOiNG DATA Af’ E PART (..IF THE PLOTTING PACKAGE 
FOR 0SL/360 S 7 M UL AT ’ ON LANGUAGE AT N P S . THEY 
• Tt - v j; r- pf rr-cmtfrP than thvs plort - ^m - g PAr»'Ar p 

7 7 AVAILABLE 



R (JT. 5 YS t N — ». 

CT GF THE VF LOCI TV VS. TINE 
PLANT ’MIS IA 


:T of TH c V.ARTATTQN CF propeller spfec 
PLAN TIM'S IA 


VS^T IMF. 


j— 


3 4.7 -220. LOO . 

IT H : PR!" 1 , PR Ml f PERC1 ,PEPC2 , PNAX 

.fi /\ N 'T 'NTS LA 


5. 


— 5 _ 


J " U . 3 4 ” 0.0 3 G . 

"T GF \ A AMO F N VS. T IMF 


5 . 


5 . 


L A l\ 1 l ' i „ S Ta 

3 l . 0.0 150. 

:r GF R,3 VS. TIME 

L 1 K N TV M V 1 K ■ > , , - - - - - 


5 . 


5 . 



"Flan i • n ■ s " i"2 

J 3.38 0.0 0.004 5. 5. 
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t:Xt'C DSL 

psl..inpu t op * 



— "S PR re- PA" ■ j A 5 rtR'TTh'i V N TSL/Vi(' S t **»Jt ATTf 
LANGUAGE . TR ANiS L AT TON CF THE PROGRAM WOULD * E 
L c OU(R c G- ■ F ANOT HE R ■ S^LLATVON -LANC i iAG F t c , t jspo 



P R OOF A A U 



LP-M i r rp 



T HF FOLLOWING PROGRAM CALCULATES THE OPEN LCOP RFSPCNSF 

■W* xppg- 5 y s T E * UNDER DYNAMIC Rtf A X TNG CCMDITTGNS. 

A LINEAR BRAKING RESISTOR WILL d E US EO , A PPRCXT M AT F r. 

FROM TUP THEORETICAL VALUES OBTAINED IN PROGRAM A . 





VARIOUS CASES CF DYNAMIC BRAKING AR 

1 1 o n ~ r> a m i>n a u c run cm i n.i ftr • 


E 


INCLUDED If: TH ~ c 














1 • — 0 Y N! A T i C B rA A K 1 1 \ G 
RESISTOR . 

R.— DYNAMIC BRAKING 


OS 1 1 iu 
USING 


U 

A 


L IN f A'K -, "V A R £ A 0 L F 
FIXED RES’ STOP 




3. -DYNAMIC BRAKING 
AIR TYPE BRAKE . 

l nv-.i am rr u a a i/ r a. n 


USING 

IICTMr 


A 

A 


MECHANICAL 

rnMK T ^AT"nN 




oe'a v ar t a si. £ r e s f stqr 

AIR BRAKE 




AND THE MCCHANTCAL 




•J.-UTUAi'ia JK 4 A IHu 

OF AN AIR BRAKE 


A l\ C A 


A 

F 


C 0 u P * K AT ?. ON 
I X ED RESISTOR . 



OFi A T THE USER HAS TO 00 JS TO SELECT A M t r H C D BY 

INDIC A TING IT BY - 1 TS NUMBER A 3 Q 

FOR EXAMPLE C AS E= 2 MEANS THAT DYNAMIC P R A K T N G 
WILL BE ACCOMPL? SHtO BY THE USE CF A FIXED 

P f ~ I T S TO R 

T N ADOTTTQN ' 



THF,, LI 

(A) 


HAS Tfl SPECIFY THE ■ F C L into 1 NG : 

- maximum permissible current :n tf^ armature 

CIRCUIT (TUTS IS ONLY LSFD '< N THE CALCU- 


(B) 


L A T T O, l r [Hu. 1 M LU^t 1 I C A L VA t u l i L» 1 Hu V ^ ^ t- . t 

RESISTOR ) 

- VALUES OP tp? FLFCTRICAL CONSTANTS OF THE 


(C> 


— NCtOR A N D GENL'TTTUR»KM"AnD KG R -s PkCT .( V EL Y . 

- VALUES OF THE MECHANICAL CONSTANTS CF NOTOP 

AMP> QC.MC 0 ATflP' , n MV. A,NP KQG R r c ° F C T T V F 1 v . 


(D) 


- VALUES OF THF INTERNAL P. p S ' S T ANC ES CF MQ-fnc 
AND G^NFRATCR, RM ANT RG 


( n ) 


—=1 — V A ! U f J S ru K — T U i! t 1 i A Lr V K Ur t L £ < * 

SF! AFT AND MOTOR SYSTEM AND VALUES FOP THE 
INERTIA OT GENERATOR RCTOR , POWER TURBINE SYSTEM 


AS I N APPEMUi < A UP 1 h t l H 1 6 1 b • 



-ETERMTNATTnH ct f PE N LOOP RESPONSE OF THE SYSTEM 

FT LG Ik a p I 
t c:lr nplot 

Jrpg r,KC 

'NST NPLOT = 1 . 
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0 E T F P " f. A T J ' 0 N CP wpr.PFLLPP PA R APET ER S 



xm rrrr i a . - v r ^ -r c.oA 



rf-T^R VNATTPM Q£ SHl£ flARAMJEJ-Eii£ 



;CKST P=2.7397E+C5 



CE TER 1 " NATION CF ROTOR CHARACTERISTICS 



"l I j S I A ^ = U . UTTTT ; H rF5 U . A 5 v 1731 XV. =TTC 0 U 5 

P E T F P M 7 \ A T T r 1 M PF GFH.= R ATiIR C H A R ACT F R 7 ST I f S 



CFST in 0 = 0. 2 5666 6 , KGG= 3 . 02 £ 669 7 , RG = C . 00C5 



t I PERT I A CF GAS-TUPfUPE/OENERATCR S FT 

<r 

TCTTSl i L 5 • 

i CTHFR CONSTANTS 



„Cf> ST PI =3. 191 59,P=2. 
f.NST OJSP-1 



initial cc net tig.vs 



INCCN iCl=30. .TC2 = 3.P3 32t ! C3 = 0 . 0 , 1 C 9 = 6 C 

4 



:tP I VAT • VF 

l-i\T r iPL ( 1 C 1 , uncr ) 



X-UTGRL ( TC3 tUJ 
t\ G T - T N T G R l. ( 1C9.GT0CT ) 



R T = T A ri F ( VI 
In =T Ar'. v ( V ) 
VF=M 



T = T A S T ( V ) 
is=yot,:*?+ (M*r )-•!=* 2 
b-SG-T I A ) 



S=N*0/B 
CT=T.\PCT( S ,VP ) 



T AsCT -"-PM r r*rt ? i -?-A 
TP= ( 1 . — T ) * T A 

TPP = (TP — RT ) / A/ 



U r CT=TPR 
V =U 

c=a 



■XX 



G G ( S ,- VP ) 



C P=CO* P * ( 0**2 )*A 
I\ = jntgrl ( TC2.ND0T) 



t t*-r 



ofH T f 



EC.: = KG-I FG l * F G T 
I A = F M/ ( P + R E> ) 



I M ; 1=(Fm- 33. ]7 (K0*P9 r ) 
Q B G =- I A * T F F * K M'* 

; 2 



P P t = ( t i\' 



4 & - R6 - 



RM 



PS R= ( ~ A**2 ) 

PRAX=EM*ta 

■ V Cl J - f ’ T-T A - I R.»" 

PFRC. i = ( PPB/P M AX ) *100. 
F E P C ? = ( PP-v /PMAX ) *100. 



T.iV = TT MF 



l-YNAP iC 
* 

r^iP.IJ. 0.2166) ITT RT h 7V 
IF(CASE.EQ.l) GO TO 100 

LE4 G A-S E. EC .2 ) GU TO 20 0 

3F(CAS c .FC.3) GO TO 300 
TF(CAS p .Fw.A) GO TO 30C 
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• F(CASF.£C.5) GG TO 300 



}nn " r- !>■'■= 4 n 



i FG=40 . 



R R=-0 .009 3 5 

1 c ( 0 tt T . 0 * 

GO TO 120 
110 RB=Q.00l5 

T2U rrrrrTTTTTH 

GO TO 60 C 



3* T 7mf+o.c17 

■ ^ rO l E ) rn . TO— VIO 



..DCCN^FO^-^A - ON ■ CBOEB - A - T-U - R -- W I O A 0 D — A — C G H 5 T A4> T PPSESrCR 
2 CO 16=0.017 

[^=40. 

i. r ; — 1 J • 

CWGT=0.0 

GO TP 600 



SE CF A'! A l R CLUTCH BRAKING OEVZCE 

j 0 U P C ^ 0 ' . T ■ ) . .2 1 2 A if ( ,M g ! 7 ) 

PS = 5.’!=TT^ 

CCL =45 3 0 . * ( R S-P C ) 

qvct = o .0 

IF(OCI. .LT.C.C) GO TO 301 

CC TO 303 

301 QCL =0 . 0 
202 CCNU'IUF. 

FV gT a l .5 0 1 : 

CC w =4550. *FH6T 

l Ft ori ,gt .on w t on Tr 

GC TO 304 
202 QCL=OCL M 

- 3C 4 ■ cr-MTT\tJF 

TF( CASE .FC .4 ) GG TO 100 
*i F ( CASF.EQ .5 ) GG TO 200 
G C r 0 6 0 C 



6 0 u ,S U M Q = 0 M 0 T - C P - Q F R - Q C L + Q M G + Q M G G 

GC TO 000 

- fOO '■F(CASG.0F.3.Ah0.CAS f i .LT.i:) GJ TO ?U 
CCL =0.0 
GO T 0 711 

TUI l>C=5.+0.2l24*(Ni**2 ) 
t< = r 'mf-6 . 7 

P-S.I=-PS-30.*K 

GCL=4550.* ( FSI-PC ) 

IF(OCL. LT.C.C) GO TO 70L 

i?fr -re 

VOL QCL =0.0 

7 J 2 C C N T " 0 1 j F 

Til QWC=O.J 

P 0 B = 0 . 0 

gRM=0. J 

PEP C 1=0.0 
?ERC2=0.0 
1 A - ■ ) . 0 
f:V = Q.O 

RB=Q.O 

>v = 40. 

* AST=-330CC . 

IW My* TFW* t AST 

„„ GC TO ROC 

3i.C SUM Q= Q A S T - C P -Q F R +Q C L 
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QFR = EC''’S'.n'<N, -22500. ,0.0, H225CC. ) 

MJU =sumo/ 12 m*2 : *1 > 

NRE AL = 6 0 . *ft 

pp o L=PRR/ l CC0C0 . 

p i ,M i «'> ! ■>!/ 1G000 « 

PVAXl=PWAX/iCCOO 0. 

IAl=* A*0.C1515 
'1J Mj i = 1 Vi /"IYC00. 

CP2=QP/10CCC. 

<- r.i i ^ iLi /i croon , 



MRL F T ft T I M= 25. ,DEL T= . 1 , DE LS = . 1 

) NT .1, EH, PRB , VOLT,OP,V, TA , PRM , PFKCl ,QFP ,ftRFAL.,P6 , PMAX, PERC? ,0'G, . . . 

ft , * AST * , oast , or tytrrstm , ■ s , l- t t, co , p s -t m. h- , tn-mtf’ , cr , p c^ -g o - F fr , - A r t , — 

I PC L , WF , ! A , OE , VP , W , T , A , B , TPR ,K , P S C 



THE c C L L f W T N G STATEMENT S ARE APPLICABLE TO tu,F 
i ~i L r ) T T * N C p AC BA QH — F - U A- PSL/3 pO . $ T - f^.'L - A - T-f Q . f — L A\f :: jAr . F 
At N.P.S . ThF ST AT FM Eft T S SHOULD BE C CINV E R T ED IF 
QTHPR TEAR THIS PLOTTING PACKAGES ARE AVAILABLE . 



CALL 0RWGUtlfTIP5,V) 
r.Ai I op not 2-lI iZILiLtJiL-AL > 
CALL DRWG< 2 PP;U ) 

‘ CALL OBWG(3,2,TIMF,PRMn 

CALL OR WG ( 3 , 3 , TI ^ , PFPC1 ) 

C ALL DP.WG ( 2 , 4 , T X M E , PFRC2 ) 
CALL DRUG ( 3 ,5,T1 M F , P»*AX1 ) 

CALL ' Oft OF- ( A , l , T* r‘ A t") — 

CALL ORWG ( 4 ,2»TT ME, EM ) 
CALL 0PWG(5,l,Tr^.DCU. ) 
Rt'.iNAL 

CALL END RW (ft PLOT ) 



CP 

C 



rrF?T7 



TABULATION OF SPEED REDUCTION COEFFICIENT LS. SPEED 



FUNCTION TABVn(V) 



T'fT TC'S VT ( i L ) , i- 1 I ( ) 



DAT A V T/ 0 . C , 2 . 5 , 5 . , 7 . 5 . 10 . , 1 2 .5 , L 5 . , L 7 . 5 , 2o . , 
hat a ,;t /£« r , r , ,!i^, . 1 1 , . 1 2 , . n?3 , . C a 5 . - 04 ~ . - 0 *3 . . ^ 0? , 



22 . 



J1£JU. 



! C . / 



GO TO l 



I P ( V . G T . 3 C . ) 

D E L V = 2 . 5 

■ ft - 1 H X ( W f t t trl — — 

SLCPW= ( to T { N + 1 ) — in T ( ft ) ) / < V T { ft + 1 )-VT ( ft ) ) 

T ABto = S LO PC z ( V-VT ( fv ) ) + WT ( ft! ) 



’FmTFI'i 

T ABC=0 . 0 

-U.TLkM 

EftC 



TABULATION OF SHIP’S RESISTANCE VS. SPEED 



-Eli NT 1 l ON T.APJUAD 

D:PFhS T Oft! V T ( 1 3 ) ,RTT(\3) 
LATA VT/O.C.2.5, 5. ,7.5,10. 



12.5, 15., 17.5, 20. ,22. 5, 25. ,27 .5, - 0 • / 



liMA a TT/ 2*0.0, 701)0. , 13000. , 2 5000. ,39000 
1136000. , 173000., 725000. ,230C0C./ 

" i m v) . g i • _ l • ) oil rn r — 

cei. v=2 . 5 

n = i e '• < i u / r £ i i 

SL0PR= (RTT (N + l. ) -RTT ( N) ) / ( VT( N + l ) -VT (N ) ) 
TAfiP=SLCPP*(V-VT(M) )+DTT(N ) 

■ PCTU^N — 

T ARP = 28000 C • 

KFTUPM 

T TO ‘ 



57000. ,00000 



i o? on c . , 



TABULATION of THRUST REDUCTION COEFFICIENT vs. speed 



FUNC T T ON T APT ( V ) 

n rm rrr r m — — 

DATA vt/o. C*. ?. 5,5. ,7. 5,10 .,12.5*15. ,17. 5*20 .,22. 5*25. *^7. 5.30./ 

04 IA T T / 3 *0-..C., .Cl , .05 . . 07 . . JP . .075 . .072 . . 3.5. .0?/ 

CEL V — 3 . 5 

N=?:ft X ( V/DELV )+l 

J. M V . 0T.30 .1-00 T r i | 

S L C P T = ( TT(N+U-TT( M ) / ( VT ( N + l >-VT <M ) 

TApT=SI.OPT* (V-VT (N) )+TT ( N) 

RE T URN — — 

TAbTaO.O 

RFTijPM 

ENT — — 



TABULATION OF THRUST COEFFICIENT VS. SECOND VC r * c IEC 
ACVANCF COEFFICIENT 



FUNC T ~ 0 N TAPC T l S . V P ) 



DIMENSION S T ( 2 .1 ) , CT 1 ( 2 1 ) ,CT?.(2i) 

DA T A ST/ — 1.,—.9, — .8, — .7, — .6,— .5, — .4,— ..p,— 
- l « 4 , » 5 , « 6 , « 7 , . . 9 , 1 . / 



OAT A CT2/-.4*-. 15, -.05, . ) 5 , . 1, .15, .2, .26, 
1. *5, .42, .35, .37, .32 , .33, .50/ 



-~ PT ' fl r, T f T- . A , - ♦ 2 . - . ? 3 , - TV?- , = TT5 — :3 H , - . -I , 
L-. 3 5, -.24, -.25, -.2, -.13, -.05 , .01, .10, .21, 

■ CFJ t S=Q 1 

N = T i- .. X ( (s+ 1 . ) /OtLS) +1 
IF(UP.LT.C.C) GO TO 2 

■ SirCTI. t (CTl (N+l ) -CT . U - M ) ) / ( ST ( N + l UST ( N) )- 



TABC t = SLCCTI*(S-ST ( N ) ) +C T 1 (N) 
GC TO 6 



■S L ' QCT 2 - ( C - T - ? ( A »l ) (ST ( ■ fV+ri - ST ( M ) ) 

TAl CT=$l_0CT2*(S-ST ( N) ) +CT2 ( N ) 

3FTIJPN 

7TT 



2 ,— . 1 , 0 . 0 , 

36, .44, .20, .21 ,.4, 

.41, — .26, — . ’’3 , — .23, 

45/ 



TABULATION OF TORQUE COEFFICIENT VS. SECOND vcr 5 T c .IED 
ADVANCE COEFFICIENT 



, FLNf.T ; flN TA8CQ ( S, VP) 

■ • r S " 0 N ST (21 » , (.01(21) ,U21L‘U 

DATA ST/-1 . , - .9 , — .8, -.7, -.6, .5, — .4, .3, .2, .1,0.,.!,. 2,.— , 

- U4 , . 5, . 6 -. . 7 ♦ . 2 , . 5 , 1 . / 

CAM CQl/-.04,-.045,-.05,-.055,-.C6,-.061,-.062,-.0fc4,-.C5P,-.04, 
i-. 017, -.055 ,-.05 ,-.0 25, -.026 ,-.015, -.003, .005 , .02, .333, .C7/_ 

D A T A '-02/- . C 8 , — • '4 35 , — • 015, 0.3, .01 , . J 2 , . C 3 4 , • 04 2 ♦ . 0 6 , • ^ 7 , .LA 5 ,.05 , 
1.062, .07 , .062 ,. 06, .053, .053 , .05, .08/ 

_ ofl S= 0 . 1 

e x t (s+i. )/Ouis) rr ~~ 



/F(VP.LT.O.O) GD TO 2 

SUlCiU ■ < CC 1 <N H J -C(, . H ( ; l ) ) /( IT ( , N4 - U - ST ( N ) L 



TABCQ=SLUC01*(S-ST( N ) ) +C01 (N ) 
GC tq ^ 
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SLT.02 = (CU2 »N+l ) -CQ2 ( U) ) / ( ST ( N+ 1 ) -ST ( (v ) ) 
T A 50 Q = S L C,G C 2 * ( S - S T ( K: ) ) <-CQ2(fy) 

i< 1 i I J K N 

LN n 



T A BULAT TUN OF PARAMETERS OF TH f P R ? M F i^PVFR 



FUNCTION T fl EO E ( V 3 * W F ) 

HiMF^s: DM V3T (5 ) . QFT ( 10) ,UET ( ICO , 100) 

V;T ( 1 J =0 » c 

V3T ( 2 I =1 OCC . 

v^t ( 2 ) =2000 . 

V ? T ( 4 ) =3 OC C # 

V3 T (5)=40CC. 

i i 3 - 3 ) C . 

h R T( 2)=43CC. 
wFT (3) = 5 3 0 C . 

!«F I ( -t ) =b3 F C, . 

5) =73CC. 

li£lXfJ-=JL3£LG-« 

.vl-T(7)=93CC. 

»>F T ( 8 ) =10300. 

WgTt-D-nsrr, 

GET ( .10) = l? ICO. 

QET U , l ) =2 CCCO . 

C c T ( 1 1 2 ) = 2 e ■ 0 C . 

Q^T ( l » 3 ) =2 7CC0 . 

O r -T ( 1 * 4) =45CC0. 

OFT (.1,3) =52500. 

QET ( 1*6) =6 2 COO . 

OFT ( ). , n^7CG0 3. 

QET (1,8) =77000. 

QET ( l » ° ) =3t500. 

o t t " ( x ; T 'TT- s -q-^oo o - : 

QFT (2,1) =12500. 

Q F T ( ? . ?) =7 000 1. 

QFT(2,3)=27CCC. 

QET (2, 4) =34500. 

QCT( - £ . 5 ) ~ 4 1 5 C 0 . 

<3ET(2f 6)=4EC00. 

QFT (2.7 ) =54GCC. 

G£ T (2*5) g 5 5 C CC. 

QFT (2*9) =64000. 

Q C T ( 2*10 ) = E S C Q 0 . 

QET ( 3 , l ) =8 CCO . 

QET ( 3, ? ) =1 TOCO. 

GE T ( 3 r 0 ) = 1 B 5C-J . 

Q ET ( 3 * 4 ) =2 4CC0 . 

OFT (3, 5) =-’3000. 

QhT(3,6 )=3 400 T)- 

QET ( 5, 7) = 3 8 c 00 . 



0HT(3. 81=45000. 
QtT( 3 * 9 1 = 5 C C C G . 



QtT ( 3, 

QET (2 , 
- 4FT (-4 t 
QET (4, 
QET ( 4, 
- QET (4y 
G F T ( 4 , 
Q E T ( 4 * 



JTTT4, 
QE T ( 4 , 
QFTU * 



10 ) =55000 
4- ) -= 5 COG . - 

2 ) =9CC0 . 

3) =12CC0. 
■T )- - ' t 1- ^00 . 
5 ) =2 2GC0 . 
6 ) = 0 6 C C 0 . 



QET(4, 
QFT (5, 



7 J = 3 COCO. 
8 ) =3 4CC0 . 
) -38500 , 



T 7ET ( 8* 
QET ( 5 , 
QET ( 5 * 



OET( 5, 
GFT (5, 
-QFT ( 5 



10 >=43000, 
l ) =3000 . 
21-65CG. - 
3 ) =10000. 
41=12000. 



C r T ( 5 , 
QFT ( 5, 



57=1 KCO.~ 
6 ) =20000. 

' ) = 2 4 C C 0 ■ 



8 ) =27COO. 
G ) =3 C5C0 . 
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ten a tio)*24ooo. - 

i F ( '/ 2 . L 7 . C . C ) f, C TO J. 

f — . i- 1 vi.bi .4urir:i gl ic 2 

X 3 - V 3 

i:r to I'lo 

X 3 = 0 . 0 
GC TQ 100 

• \2 - 4pOO. — — — 

GC TO 100 

IOC IF(Lir .LT.3 3CO. ) GO TO 3 

. ftftftti rrerurr: i uj m ^ 

WF1=WF 

c-r to ?no 

1 W P l =33 00 . 

GC TO 200 

, n F l - - -F 2 3 - TO . 

2 oo i=:r:x( x3/iooo. )+i 

J = I F I X { (WP 1-2300. 1/1000. )+l 

3 k -X 3 - V i T (i 1 “ “ 

CVs^FJ.-WFT ( J) 

LLEi .o; = (np/inno. )*<qct<’ =h ..n-urn T,.in 

CtLQW=(OW/lCOO. )--MOFT{T , j+d-ofT ( {,j) ) 

TA6GF=0ET( T , J ) +C ELOR+OELQW 

GFTU P N 

EKD 

rcrr rr tlC ' A T NG “WTr TVT F P A F~ T L T V — f^F ' P rtTTT * "! 3 P AClwV-e 

c 0‘} DSL/TfcO STMULATTON l ANC-UAGF AT KPS. THEY 

SHOP! H P F f-'O f T F T F 0 7F OTHER THAN THIS PLCT T T|.G PAC.<AG C 

T S A VAT LABI. E 



'PLCT.SYS T N 00 * 

.GT CF THF VELCCTTY VS.TIf‘5 

.FLAMIN' *, TA 

.0 6 . - 2 . 7 . 5 . 

rj f F THF VAPTATTQM or- OPOPFI I FR SPFFP VF.TTMF 



FLAMIMS 14 

C 4.7 -230. IOC. 

CT CF PPiU , Pavi , PE H C 1 » P FP C2 » P MA X 



FLANT7MS 1A 

0.74 C.O 

CT CF rf.- T VKC 1 E K"VS ' . ' T 1 ' f- ’ P 



8 C . 



FLAMIMS 1A 

I tf SB vitTTTE 



C. C 



150. 



FLAMIMS 
■Q 0 . 38 



1A 



-O^G- 



0 . 00 4 
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FXFC 0 r > 1. 

■ .input jd * 



— Eii£H£AM HAS hiLULLEN i ll USLZ1AQ SJJZLLULLLC& 

LANGUAGE . TRANSLATION of the PROGRAM WOULD BE 
REDO! RFC T F ANOTHER SIMULATION LANGUAGE T S USFD 



TTc FT r nT R T TT? — r 



ld£ r n t rviruo prugp.au caicuiatfs thf npru ippp 

n c t«e system undfr dynamic braking conditions. 

THE INERTIA OF MOVING PARTS IS BEING USED FCR 
DYNAMIC BRAKING. — 



B-E-S£HLlS-E 



TU r US 


C P HAS TO SP C C T FY ThF FOl FL'mTNC : 


( A ) 

( 


- VALUES OP THE ELFCTRICAL CONSTANTS CF TH r 
MOTOR AND GENERATOR, KM AND KG RESPECTIVELY. 

— -VAIIIF'" f ,r: tl.t Mcru>l:1f ai rnMCT aAiTC rc ur. me 


( C 1 


AND GFNFPATOR, KMM AND KGG RESPECT ! VELY . 

- VALUES OF THF INTERNAL RESISTANCES OF MOTOR 


( D) 


A i\ fl u E N E K A T G t R m A N C Kb 

- VALUES FOP THE INERTIA CF PRCPELLER , 

SHAFT AND MOTOR SYSTFM AND VALUES FCR Tmf 


INERT! A OT GEN ERATLK RCTCR .POWER TURBINE SYSTEM 
AS IN APPENDIX A OF THE THESIS. 



CP.S T? =C t I N * T I A L 'T S H S T HE CRASH-STOP COMMAND 

Cf.STP = l , PRCHIRITES THE APPLICATION CF THE C MEPGENCY 

PROCFOllRFS ANI THE SUP 0 P E R A T F S AT CCN P "7 T CNG 

REGULATED BY THE INITIAL CONDITIONS AMD 
STEADY STATE DYNAMICS. 

TEG TRAPZ 
TGER N P L 0 T 

T.GFk OPSTP 

NST NH.UT=l 

G^T-g.PMTNATinN OF PROPEL Lb-R PARAMETERS 

NST 1=30171. ♦ C = 1 3 .9 4 

CETFRMIL ATTON OF SHIP PARAMETERS 

TTT' M = 2. i'3s7l- + G5 

t-Ii-I-S-B-N T A A T I O \ OF MOTOR. CH A R A C T ER T-S-T. I CS 

\ ST K,M=Q. 06 3 26 t KMv=Q .459153, RM = C. CC0 5 

DETERMINATION OF GENERATOR CHARACTERISTICS 
TTT K C=D . 7 b t Z S T, 7 K'GG= 0 . U 2 * A 6 A I, k < i“c TTHT C 5 

— Tf.TQTJA n F iTAC-TKRnrMr/^r-JFftATf a S F T 

V| ST I C- = ^ C1 1 S . 

C T HER CONSTANTS 

I**! P.i= i. I'43.b9 ,P-2. 

iA-IT'AL -CC N C I XI PALS — 

,-CN ICi=30.,IC2=3.6333, T C3=0.0 ,'C4=6C. 
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ST CRSTP=1 



U= lM1Rt. ( TCI * UDC T ) 

>=• NTGRL ( T C2 f U ) 

( VT 

k= T ABW(V) 

VP=( l.-W )*'/ 

— i = • rrrrrnn 

a=vp** 2+ (r;«n )**2 

P = r ' W r~ UJL1 

$=K*D/B 

CT= TABC T ( S »VP ) 

T fl - 1 ■< ■■ n ( C ' ■ > - i? i v A 

T P = (J . - T ) * T A 
TPR = ( TP-RT ) /M 
urr.y= rpp ~ 

V =u 

C C = T / li C (J ( S »VP ) 

.CP=CC*P* ( r s> * ? ) * a 
NG 7 ='NT 0 RL ( TC 4 .GTD 0 T ) 

N - T N T 0 a L (TOP i NHCiT ) 

1FM=40. 

I FG 1 = ( C w -3 3 . )/( KG*NGT ) 

CVbn. r ji - - f Fvg iTF gAL 

EG=KG*T FGl*NGT 
lArlf-.-FN 1 ) / ( EG ±£±1 ) 



44-0 



!F (CRSTP. EC.O. ) GG TD 110 
1 F ( Cf S T P . L C . 1 . A NIT . f ' ! « t T « 0 . 2 1 MS ) GC TO ttn? 

Vm F — 2 GO 0 « 

GC TC POO 

GO WF = 3? )0.+lCCC.*RA M P (1.4) 

I F .< ur ...C E -. 12200 . ) CO TG 1 10 

GC TO 1PC 
1C V«F= 1.230C . 

Z C CCMT'MIir 

CC C r, M T T NIlF 

Q V = K*>'H*7 A * T F M 

CGF*'-KGG*t A*lFGl 
GF = TABO c (NG »WF ) 

— n ;F-r:c-Qr,FA 

GTCCT=D!F/(?.*PJ*IG> 

N3=NGT*60 . 

SUEC-OM-CP-C F R 

CFP = FCiVSMG,-22 500. ,0.0,+2 25 0C. ) 

NDCT = SIJMC/ (?.*P T*I ) 

NRfcAL = 60.*K 



PL E 

TRL FT Ml M=?5. .CELT=0.001 ,DELS= 0.0125 

If f . 1. 1 3 t V t Q P t ' ! F 0 1 1 M G T « C T > P T rfe^Hri A . N W » QGCN t QF , Xf - 1 A , T , T PR , ♦ . . 
Cl F. f EM , 6 » TA , TP , N3 » EG * SUMO t GF P , NDCT ,NR E AL f X » I FG2 . I FG t Q C » P F . EF 



-T-H T . ■ FO L LCVIMG STATEN EM T S A P F APPL t CABLE TO THE 

PLGTT TNG PACKAGE FOP 051/360 SI M ULAT TC'N LANGUAGE 
A T m.P.S . THE STATFMENTS SHOULD BF CONVERTED If 
OTHER t h an TU i- S— P tr OTT 1 ! NG PACKAGES ' AR ~€~ "A VAT L Apt. F . 



HCI 0 R L G ( 1,1. ,TiMb,V) 

CALL 0 P W G ( 2 » 1 ,TIPE,NREAL ) 
-C-Al L n R W C { 2 » 1 »X J P E » T F G 1, ) - 
CAL L DRWG( 2 »2»TTME» ! FG ) 
CALL DRWG(4»1»TIME»N3) 



128 





CALL 

CALL 


OPWG( 5 ,1 ,TI ME, RC) 
dp w r - ( 5,2,t;ne,rm 




mnai 

ft! 1 


p Nn runpi nr ) 




p 




• 



TRAN 



• TA8UL A - T - IPU - OF S,[?FEP REDUCTION C Q c F F V C T J-- VS . ^ R P 3 D 



Iwrtrrttrrc rmtr> 

Dlf'FNSI CN VT (1 3 ) ,VvT (13 ) 

liAIA VT/O.C .2.5 . 5. ,^.5 . 10. .12.5, 15. , L7. 5 , 20., 22.5,25. , 27.5,30. / 
U AT A WT/4*C.C,.O05,.Ol,.O2,.C23,.C45,.CAf>,.03t;,.02,.0CA/ 
l F( V.r,T .30 . ) GO TC l 

Oil w-g . 5 

K=?P iX(V/CELV)+l 

Sl.CPW=UT(N+1. )-wT(N> )/( V T ( N+l ) - V T ( N ) ) 

Ttrrr?r= SL^- P W ‘ (V-VT ( r, ) ) mr (m ) - — 

RETURN 

TABH-Q.Q 

RFTURA 

END 



TABULA TU T" ' "O* 1 SHIP'S RESISTANCE V S T^P^t 



FUNCT'i ON T A B R ( V ) 

BIVENS’ ON V T (13) » RTT ( 13 ) 

DATA . y-W 0 . C , ? . 5 t 5 . , 7 . 5 , TO . . l2 ,. 5 -«- l 5 . , l 7 . 5 . 
DATA RTT / 2 ^ C . 0 » 7000. , 13000 . , 25030 . ,39000 
1136000 . , 172000. .225000. ,280000./ 

U ; ( V.OT .30 . ) GO TO 1 1 

C F.L V= 2 • 5 

K=J F iX ( V/DELV ) + l 

-TTUPI-— TR T T ( N + I ) -P T f { T; ) ) / ( V r ( N+ 1 ) - V T ( v ) ) ■ 

T AB R= SI.F) PR * (V-VT ( N > ) +RTT ( N ) 

- RFTliftM 

TAP.R = 23090 0. 

RETURN 

— trfr© — 



7 0 . . ? ? .5,?5. ,2 7.5, 30. / 
,57000. .SCOOC. , 103C00. , 



TABULATION OF THRUST REDUCTION COEFFICIENT VS.SPFET 



FUNCTION T ART (V ) 

"nil V'vt/q.cIoH 10. ,12.5,15. ,17. 5, 20. ,22. 5, 25. ,37.5,30./ 

DATA TT/5*C. .C..C1..05, .07, . 03, .075, .072, .06, .02/ 

~CbLV= 2 . 5 



N = I F f x ( V / 0 E L V ) + 1 
AF( V.QT.3D. ) GO— 



XU— L 



SLOPT= ( TT( N+l )-TT ( N ) ) / ( VT( N+U-VT (N) ) 
TABT=SLOPT* (V-VT ( N ) )*TT(N ) 



TABT=0 .0 
RETURN 



TABULATED OF THRUST COEFFICIENT VS. SECOND PCDIFIEC 
ADVANCE COEFFICIENT 
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rrnnTTnr t jftth s , w i 

r a,i r m c r n k: ct / 01 i r t i 



DATA CTl/-.4,-.2*-.23»-.32,-.35*-, 
1-. 3 5, -.34, -.25, -.2, -.13, -.0 5, .01, 

C.U 5 = 3:i 

N = IFTX( (S+l.) / D G L S ) + 1 
•'F(VP.LT.C.O) PC TO 2 


.38, .4,— .41,— .36,— .2 C ,-. P 5, 
.10, .21, .45/ 


SI.0C U= (CT KN + .U-CTl < N))/(ST(K + 1.)- 
TABCT=SLCC U*(S-ST(N) )+CTl (Ni 
o c T 0 ^ - - 


- S T ( N ) ) 


SLOC.T 2 = ( C T 2 ( N+l ) -C T2 (N ) )/( ST IN + 1 1 - 
' TABCT=SLCCT?4(S-ST(N))+CT2(N) 

i RETURN 

EMC 


- ST ( N ) ) 


TABULATION OF TORQUE 
ADVANCE CO EFFICIENT — 


COEFFICIENT VS.SECONC MCCI c TEC 


KTj'JUJf'r TAHCu(StVP) 

0 I KINSTON ST(21),CQl(21)tC02(2U 
G A T A S T / - 1 6 i * • 5 t ~ < 

i • 4 t • 3 « .6 t *7 t .Bt »9 t l • / 


• 4 ? ~ t“»2 *“• I. ♦ 0 • » • 1. 1 • c * • ? » 



DATA C0 1. /- .04 ,-.045 ,-.05,-.055,-.06,-.06l,-.C62,-.C64,-.C58,-.04, 
■ m? . .. .. ..as - 03 s . a . - . a i s . - . nns . . ans ir ? m / 



OAT 4 CQ?/-. 08, -.035, -.015,0.0,. 01, .02, .03 4, .04 2, .06, .07. .045, .05, 
l.C6B,.')7t.C6i,.Ofct.053,.05 3».C5».Cfc/ 



ur.L S-'J • 1 

N = I F I X ( ( S+l. )/DELS) +1 
•:f(vp,! t.o.o) on tn ? 




SL(J0Q1 = ( CO l ( N+l ) -CQ1 (N ) ) / ( ST ( N + U - 
T ABCQ - SLOC C 1* ( S -ST ( M ) ) +CQ1 (N) 

C-F TO 6 


-ST(N).) 


SL CCQ2= ( CQ ? (N+l ) -CQ2 ( N ) ) /( STIN+1)- 
T A8CQ=S LCC 02* (S-ST (N ) ) +CQ2 (N ) 


-ST(M ) ) 


■ KT.'T U ,(> N 

ENO 





TA BULAT ! ON OF PARAMETERS OF THF CRIME MOV F R 



FUNCTION TARQF(V3,WF) 

DIMENSION V 3T ( 5 ) t WFT ( 1 0 ) ,QET( 100,100) 



— vjim-u.u 

V3T( 2) =1CCC. 

V3T (3) =2000. 

V 2 r ( 4 ) = 3 0 0 C . 

\/3T(5)=40CC. 

WFT ( l > =3 30C 

W FT ( 2 ) =43 C C . 
wFT (3)=530C. 

fcTTK) -63CC. 

W C T ( 5) =7300 . 

WFT ( 6 ) =8300 . 

(7) =<530C. 

WFT( 8 )=102C0. 

mJT(Q) =U3C0 

WFT ( 10 ) = 1 2 3C0 . 

QET(1,1)=2CCC0. 

■ eirT (1 ,2 ) - 26500. — 

OFT ( ]. , 3 ) =3 7CC0. 

GET ( 1 , 4 ) = 4 5000. 

Ohl l 1, 5)=535CU. 

0ET( 1,6)=62CC0. 

■ 0£x n. ,7i = 7croo. 

0FT(\ t 8)=77CC0. 

OET (1,9 ) =86500. 


r 
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OPT ( l. 10 ) = <55CCO. 

<;!:!(?, 1) =12500. 



I j>- ' U. 1 T*7 er r 
0 5 T ( 2. , 3 ) — 2 7 C C 0 . 
opt L2-tAi=3-^££ 



U r I 1 £. * » > ~ V.: u * 

0feT(2t S)=A1SCC. 
CET(2,6)=AE000. 
r. i PT l ? ■ y - > -5/ t OCC. 



QFT (2» 3 )=5SCC0. 
QFT(2* r >)=6ACCU. 



UFTty.t 10 ) -t * 000, 

0 FI ( 3.1 )=. q con. 

0E112 ,2) -130C0. 



(OFT (3, 5 )=l?GC0, 
G F ; T ( 3 * A ) = 2 A C C 0 , 



i a fT ( 






UFII 3, 6)=3A0CC. 

Gf T (3,7)=3S5C0« 

crnorsi-^cco. — 

OF T ( 3 * 0 ) ■= 5C00C . 
(OFT (3.1 0 ) = 5 5 C 0 0 . 



OFT ( A, J. ) = 5CCC . 
(OPT ( 4, 2 ) =<?f CO. 

o r- t < a , 3 ) - 1 ; c o (/ 



CFT(4,4)=17 C 00. 
OFT (A, 5 ) = 2 2 C C 0 . 



(OF T (A,7)=?CroO. 



=:-acou. 



0 1: T 



OFT (A.. 10) =43000. 
OFT ( A , n = 3 000 . 



0 FT (£■»?.)= 6 5 C C 
CFT (5,‘=‘ ) =ICC00. 



oft< a, A )=i^ - e -e o -r 

CFT ( 5 » 5 ) = 1 7C00. 
OFT (3,3) = 2 C C CO . 



QfcT(P, 7 ) = 2 A C 0 0 . 
CFT (5*8) “27C00. 
0 FT (5 , A) = 34 SUM. 



GET (5flG)=?AC00. 

IF( V3.LT.C.C) r -C TC 1 
-fPt 3 r^T^ T 7 4CCC. ) ' GO TC ? 
X ? = V - • 

c.r: tq too 



X3 = 0.'J 
GC TCI 100 
-x.3=Aa.%o, 



GC TO IOC 
ICC IFIUF.LT .32CO. ) GC. TO 3 



.1 ' ( <•' i~Tj 

Vi P 1 = V F 
CC TO 200 



GT .12300. ) GO TO ' #- 



'vsFl = ’300. 
GC TO 2 OC 
-l»Flsl2 300. 



CO I = 1 F 7 X ( X3/ ICCO. ) +1 

J = T F T x ( ( W F 1 - : 3 0 0 . ) / 1000. )+l 
— eff 



X3 - V3T ( T ) 

Cu = un.-WFT ( J ) 

DFLCF'=(DP/\OQO. )*(QFT( T+l * J ) - U E T ( 



t J) ) 



TELUU=UJW/ l COG. i v IGF I U , J + i J =G U l 07071 
T A F Q E = Q E T ( I , J ) + C Et.QR +DELOW 
-F.CTLIF V 



tNC 



- T - l 'f C r n L L C W I H G ' D ft T A A P-F 
FOR 0SL/36C SIMULATION 
SEOUL T EE MOCTFIED T F 
TS A V A I L * F L E 



■ PART or T I C PLOTTING ■ P - ACK ^ fr F 

LANGUAGE AT N P S . THEY 
C TF ER THAN THIS PLC T T T NG PACKAGE 



U r T . c; V £ r M p n ^ 

■ T CF THE VELOCITY VS. TIME 
‘LANT.T N T S 1 A 
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?T LF T FF VAPTATU'M OF PRO°FLL ER SPEEC VS’.TIPF 

fi at. i :n'"s ns — 

c 4.7 -230. ICC. 5. 

-t l£ p d -I . p? m . PFp.n . pr?r.? .p vax 

FLAMIN' S 1A 

I) 0.34 0.0 80. 5. 

-4— -t-f ' A 4AFP — R* VA.T ;wr 

HAMTNIS IA 

0 l. O.G 150. 5. 

C l cr R . R- -VS. TT .v r 

FLANTtNTS 1A 

C_ Q . 3 5 Q^lI} 0.004 5. 



5 . 



5. 



5." 
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// EXEC DEI. 

//tSL. INPUT DO * 

T 



WS orcgra.-waS WTthn in QS)/R h o G~ MUI AT TIN. 



, ANGUAG E . TRANSLATION C F THE PRGGRA V WOULD - BE 
REOU T RED IF ANOTHER SIMULATION LANGUAGE 7 S USED 



(i 



* 



rmr 


£ 


* 


THF FOLLCHTNG PROGRAM CALCULATES TH- Cl OSEO LOOP RFSPCNSF 


r 

* 


OF THE SYSTEM UNDER DYNAMIC BRAKING CONDITIONS. 

WHEN THE INERTIA OF MOVING PARTTS OF THE GENERATOR 

AJ1LO X Jrj nn jii: q TJifl.ftl MP >\ J £ ii_cg_n QJ2Z LiVM.ni T r gj ait T Mr 






r~ 

i' 


the IJS^P HAS TO SPECIFY THE FOLLOWING : 


(A) - VALUES UH THE ELECTRICAL CONSTANTS THE 

:< ‘ / f, TOP AND G p NF,< A TOR , KM AN C KG RESPECTIVELY. 

1 l t 1 = \< a 1 1 1 r- <: n p n-u: m f r m a v t r a , r nv <;t a r-T<; re v . 0 r ~t p 


't 

f 


AND GENERATOR ,KMM ANIC KCG 'RESPECTIVELY" . ’ 
(C) - VALUES OP THE INTERNAL RESISTANCES CF MOTOR 


\ 

t 


A !\ 1. GL Ml K 4 ] IJK , R ;•! — A X L — rV. 

(0) - VALUES FOR Thp INERTIA CF PROPELLER , 

SHAFT AMO MOTOR SYSTEM AND VALUES FOR THE 


Jl 

•( 


INERTIA U! ( 'L N ERA T LR K L TOR , RlJwL R TURBINE S' Y S T E l/ ' 
AS IN APPENDIX A OF TEE THESIS. 


’< 


CP S T p = C , TNI TI AL' f S C S THE CRASH-STOP COMMAND 

CRSTP =1 , PROH T BI TES THE APPLICATION O c TH E EMERGENCY 


* 

(• 


Jufnrj^ c Is AM', ( rc UKLKAIro A I UJ"J *! !<_• o 

REGULATED 1Y THF INITIAL CONDITTGNS AND 
STFADY STATE DYNAMICS. 


> 


(NT EG 
(NIGER 


T'PAPZ 

NPLOT 


N T : j ir u 
: r j(\ST 


CRSTP : 

NPLOT = \ 


s 

X&ST 


n.- T HRiM 1 NAHO'N""CT- PROP ELLER PARAMETERS 


-W.Di.t7J. . . r - 1 7 . G 4 



GET r -R M ! N A Ti ON OF SHIP PARAMETERS 



IONS!' M = ? . 7 3 '5 7 E t C 5 

I'Ll r: P M ", T\ A | A I J N OF Ml!! I OR CHAR AC I EK 1 S I • CS 

XkSI * v = 0 - IVR V7A t If Mjj= - r i '- ! I •> ~ . 1M=TU-Q£L15 

OETE RM ’NATION OF GENERATOR CHARACTERISTICS 



ON ST KG=3. 2 36 66 6 , KGG= 3 . 0 7 36697 , RG=C . 00 C5 

T T^FRTTS ET GAS- I UL h i NE/ GLNLh A Tuk STO 

7 

XA-S-I KU33X6 

t 

f OTHER CONSTANTS 



IfNST PT=3.1415S ,F = 2. 

} 

; O'l i 1 7. A L c L A C T TTU N S 

XT pn I r 1 =3 "I . . r r, 2=3^ a a a 3 . ■ r. i=g . 0 . 7 r. 4 = 6 0 , 
Xi\ST CKSTP-l 
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I-RAT MJMf l J * 0 fc N ( ? ) » C l ( l J » f 2 ( L J 

r.L F M |M ( l ) = IP . 7fcF+ Q4 » DF'I( 1-3) = 1 . ,300. ,19. 76F+04, C 1(1 ) =10 . ,C 2 ( l ) = 40 . 



pyy/aT n/c 

X = " NTGF L ( T C3 ♦U J 

__ [ T-TAkk ( y ) 

l\ = T A B W ( '' ) ' 

VP=< l.-'-J > V 

t = t rtt i 

A=VP**2+ ( N *?. ) * * 2 

-< = S0i T ( A) 

■ S=N*l)/ r B 

CT= T Aif.T( S ,VP) 

TA-GT .t p a <C**2 > *A 

TP=(1 .-T )*TA 
TPR = ( Tp-RT ) /vi 

urrT= Trr . 

v=u 

r:c = T A^r.o( r » v p ) 

CP=C(>*P* ( 0*?)*A 
NGT=IMTGPL ( T C4 * GTDOT ) 

G-vntgrl ( tg? PW ) 

:fm=4o. 

\r-C,-[ . ) / ( kg*ngt ) 

r : i- = h-i:^is'.-' ( i fg — » - h j o . , u . c , f 4 o o . ) 

kC=FF/ T FG 

k F = T R. "1 F R (0 . ,2. , C I. ,MJM,nEM , RC ) 

i c 0Ut-/H 

rP = k^*I F m* k R F AL 

If c = K. M v F- r. 1 ^ A -G T 

!A = ( c G-F« ) / ( PG + RM ) 



4 AM 1C 

irtCP'fTTPTPC.O.) GC T C 110 

rF(CRSTP.bC.l. AND. fJ.LT. 0.2166) GG TG 10G 



WF = 3?,00. 

GG TG ,200 

LOO 4 1 ' - 3 J? 0 0 . + 1 C C 0 . * is A m P (1.4) 

£F(W C .GE .12300. ) GG TG HO 

GO TO IRQ 

LlO vU-=t2?00. 

120 COM T r NU e 

2-00 f.Ck: T? M'JF 

GM = KMM*{ A4 TFM 

0GEN = KC-G*t A4TPG1 

ivF^TAOQL i k ' 2 - ^ fi 

CtF=0E-0GEN 
G t n 0 T = 0 T F / ( 2 . * P 1 * I G ) 

I'M J=NGI *6 0. 

SO^O-O'-i-GP-CFR 

U££-=JLCiiSliI iv T -22FGn. f P-H. +?^SHO. ) 

i\ 0 C T = S U > M Q / <2.*PI“I ) 

K4FAL=60.*k 



FTTn * W=Z5 . , U t L l = J « ) J r t iJET 5=U 70 1 4 P 
. I , S, V, QP , TFGi ,NGT,C T , RT , GM, J A,N , C 
n - e . ~ a - r. . t a ■ tp . lb . c a . sim . qfr . NH L 



CQ 



. nr.f;j . op 



T! i r F 03 < LC ' »l T NG ST ATF - HFNT G ARC A PiH ,'tC A PL F TG T‘ +f 

PLG T T ; \ G PACKAGE FGR JSL/3fcO SIMULATION LANGUAGE 
AT \|.P.S . THE STATEMENTS SEOUL D BE CPNVFRTEO IF 
TTFTFP — IT-AN — TH L S PLUMING PflLKatbS AKF AVAILABLE 7 



CALL QRWG( 1 . f l,TI*F.V) 

CALL DR WG ( 2 , .1 .T I y E, NREAL ) 
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) 







CALL 
CAL L 


0°WG( 3 , l , T I vi p , t FG1. ) 
D R W G ( 3 ,2, TIME, I FC5 ) 


1 p Z 

CALL 
r ai i 


DRUG (5,1, TIM* ,RC ) 
OR HI ( 5 .P.Tl MF.RF) 


wINAL 
L *(A- 


t'.NNRV, ( N PLOT ) 



p 



Tk A i\' 



TAi-iULAT * OM O F S ' P E i: C HC & d ' CT f QM COEFFICIENT VS^PrEU 



FUHUTTrTTJ TABMV) “ 

DIMENSION VT ( 13) ,niT ( 13) 

Q vr A V T / D ■ C T ? , 3 , E . 7 S . 1 A . I ? , § . 1 5 , T 17 . p T 2 '-) . , ? ? ■ 5 ■ 33 . , ?7 . 5 , 30 . / 
DATA JT/4*C.C*.005».01».02,.023,.Ch5».045,.038*.G2,.OC4/ 

J F( V.fiT .30 . ) GO TO l 

ep.v-r . r j — — 

N=! -IX( V/DCLV ) +1 

SI. CPN= ( WT ( N + l )- 4T< N!) ) /(V T (N + l )-VT (N ) ) 

T\BVi=SLOI i W!' (V-VT ( N ) ) + WT (N) 

RETURN 

TArU=/U3 

RETURN 
t n n 



TABULATION O'- Ship* < RESISTANCE VS. SPEED 



EbNOT' ON TAfiR(V) 
DIMENSION V T (13) ,RTT( 13) 



" U A T A 7 T/C.C 



DATA RTT/2*C.O, 70 00 . , l 30 00 . ,2 50CC . ,3900 0. , 57000. , 8000 C . , 1C3C0C. , 
1136 )03 . , 173000. ,225003. ,23 0000./ 



‘ 7 - . - 5- , 30./ 



I 775 I .30 . ) GO 
C t L 7 = 2 .5 
rTc • v r \i / n i 



m i 



-U 



X ( VvLQ-EUg-U. 



S10°R=(RTT (N + l) -RTT (N) ) / ( V T ( N + 1 ) - V T ( N ) ) 
TA6R-SL0PR* (V-VT ( U ) ) + RTT ( N ) 

-h CT I JRM 

TARR=23000C. 

RETURN 

-pjTl 



TABULATION OF - THRUST R EQUATION COEFFICIENT VS. SPEED 



FLMCTT'IN TAPT(V) 

■ ' PTETTS FJN VTTT3T7TTTT31 

DATA VT/0.0,2.5,5.,7.5, 10. ,12. 5, 15. ,17. 5, 20. ,22. 5, 25. ,*7.5, 30./ 

oat a TT/^.->r f n,.ni, >H5 ♦ -n 7 , . g a t . n 7 s . - 07? .. On f . fl ?/ 

OF.L V =? . 5 

N= I f y X ( V / 0 E l. V ) + 1 

tH 7 . GT . 30 . ) G 0 -TO" t 

SLCRT= (TT( N+l )-TT ( N ) )/ ( VT ( N+l )-VT (N ) ) 

TAB T=SLOPT* (V-VT (N) i +TT ( N ) 

ITT- 1 IJRN 

TA6T=0. 0 

-RFTIIPM 

ENO 
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TABULATION OF THRUST COtFFICTFUT VS.SFCrRT FTr^F TC r 
AOVANC E IXEFF I U.LNI ' - 



FUNCTION t APCT(S,VP) 

D(f-* c NSTnN ST (21 ) , CT1 ( 21 ) ,CT2< 21 ) 

— — 1 — •> T / . — A. » . — « ' J i * i '! i - t7 , . o , - . 5 , — » 4 i » ■ 3 . 2 

1.4, .5,. 6 ,.7, .8, .0,1./ 

HA ? A CT2/-.4,- .15,-.05 , .Ob , . I , . 15 ,.2, . 26, .36, . 
"T • I i • T 2 t • 3 L r .37 * • j 2 » • _* 3 * • L / " ~ 

CAT A CT1/-. A, -.2, -.2 6, -.32 ,-.35, -.38, -.A, -.41, 
l — . 34, — .2^, — .2, — .18 , — . o 5 ♦ . 0 i , .10. .21. .45/ 

0 h L 5 = 0.1 

N = I F T X ( ( S + l . ) /OF LS ) +1 

,:r-< JU . LT . O . 0 ) fi- 0 TU 2 

St CCTIMCT I (N+1 J-CTl (N) I / ( ST(N + l )-ST(N ) ) 
TABCT=SLCCT1*(S-ST(M) )+CTl ( M I 

sfi— m n 

SLCCT2= (CT2 (N+l ) -CT 2 ( N ) ) / ( S T ( N + l J-ST(M) ) 
TARCT-S'.nr. 7 2* ( S-ST(N) ) +CT2 IN ) 

RETURN 

END 



4 ,- C . 1 , . 

44, .3 3, 

- . 36 , -. 



i t » /. t * -/ t 

. - i , . 4 , 



TABULATION ()r TOKOU E CCEFFTC LENT VS. SECOND RCO^IEO 

rtOVANC i : CO EH iCIfM 



FUTT7T! iUN i TTFTCTJT S , VP ) 

C I MENS T ON S T ( 2 1 ) , CQ 1. ( 2 1 ) , CC 2 ( 21 J 

DATA .ST/-1 . .-.9 , - . 3 7 .-. 6 5., -.4. -.3, -.2, -. 1 ,0.,.!, 

1 . 4 , . 5 , • o ♦ .7 , • 8 , • F , i • / 

DATA CQ1/-.C4,— .045,— .05,-.055,-.C6,— .061,— .062,— .064,— .C5F,-.04, 

1 - .IT37 7—rd*r* , - .00 , — ♦ 0 ? 5 , - . 0 2 o , — : 0 l j , - . 0 0 3 , ♦ 0 OB ♦ - .02 , . 0 3 ? , . C 7 / 

DA-M C02/-.0A ,-. 035,-. 015 , 0.0, .01 , .02 , . 034 , . 04 2 , . 06 , . 0 7 , .04 5 , . 0 5 , 
l.G63..07..C62..Cfc,.Q53,.055,.U5,.C6/ 

DFLS=0 . 1 

i\ = tc? x( ( s + l. ) /DELS) +1 

yp j. t.-Q.c ). .r,q to ? 

S t. CC 0 1 = ( C C 1 IF + 1 ) -C01 (N) ) / ( ST (,\+l )-ST( N) ) 

T APCQ=SLOCG L*< S-ST( N ) ) +CQ1 (N ) 

"SC - T 0 ■ 6 

SL 0C02 = ( CQ 2 ( N +1 ) — CU2 ( M ) ) /( ST(N+ l)-ST(N) ) 

TABC0=SLCCC2*(S-ST(N) J+CQ2 (N) 

RE I IJRN 

ELD 



TABULATION OF PARAMETERS OF THF PRIME M OVFP 



FUNCTION TAPCE(V3,WF) 

DIMENSION V?T(5) , Is F T ( 1 0 ) ,QET( ICC, 100) 
V 3 I ( I 1=0.0 
V3TI2) =10CC. 

V -T( 3 ) =20CC . 

V ? T ( 4 ) =300 C . 

V3T(5) =4Q0C. 

-x T-T ' t 1) =~:cc. 

WFTI 2)=430C. 

RFT ( 3 ) = 52 C C • 

WEI l 4 J = 6 5 U L' . 

WFT( 5)=73CC. 

h ft ( f i = r ^ r, r . 

WFT (7 ) =93CC . 

W p T(3)=lC3CC. 

tfilv) =11300 . 

WFT(10)=1220C. 

uet(i, i,) = ?cc:o. 

TJTrTTTTT^TETCTFT 
Q FT ( 1 , 3 ) = 3 7 C C 0 . 

!3 t: TU .4 1=45000. 

QFTU, 5 1=52500. 

OFT ( 1,6) =62000. 
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0!~T(], f 7)=7CCC«j. 
GET ( l , r , ) =7 7CP0 . 



m-rr 



l = F t. *. C ' > 



ivi-ut: 

QET( l, 10 ) =3 5000. 
LI £ T (P.1 ) = 1 ^03. 



0F T ( ? » ? > =?CCC0. 
Q f- T ( 2* 2 ) =27CC0. 



t * 'M ; * t jji 

QFT(2,5)=415C:. 
UPT (2 t G)=4FCC0. 



) - 3 4 m:o . 



CfT ( 7. , T > = *4000, 

(OFT ( 2 * 6)=5SCCC. 

till ( c > ) = 64000 . 



UET(2.lC)=c3COO, 
CFT(3, U=SCCC. 
< ;g icon . 



QFT(3,3)=IF5CC. 
GET (3. 4)=?4CC0. 
G r T ( ’■ » 5 ) = 2 7 c Z 0 . 
QET<3»M=?4CC0. 



J p T ( 3 t d ) = 4 500(Jv 
QET(3,9)-5CCnO. 
OF ■ ( 3 . 1.0 > -55000 , 



GFT( 4, 1. ) — 5CCC. 
QFT<4t?)=<5CCC. 
GF T < 4, 3 ) iil?gOO. 
QET(4, 4)=175CC. 
CFT(4.5)=~?P0( 



GET ( 4 f 6 I =2 tcc2. 
GET (4, 7)=2CCC0. 
GT- - T ( 4 - r £- > -=- ? - 4 COO . 



QFT(4,3)=365GC. 
CFT(4, 10)=43000 
GET 0,1 )-3CC0. 
OFT ( 5 f 2 ) = 6 5CC . 

Q F T ( 5,3 =ICC00. 



OFT ( 5 ,4 ) =13000, 
CET(5,5 ) = 1 7COO, 
0 F T ( S f fj ) = 7 C 0 0 0 , 



U r T ( 5 . 7 ) = 2 4 C 0 0 . 
OFT ( 5, 0 )=27CCO. 
Gl ' T 1 r ' i 9 )- - ? e CtF r1- ;— 
OFT ( 5 * 1 0 ) = 34 COO . 

T F ( V 3 . L T .0«0) GO 



(V‘ .GT.4C0G. ) 

X 3 =V 3 

-G.T — T.Q — 1-0 0 



TO 1 



TUTTUr 



X3=0.0 
GC IP TOO 
"X 3 - 400 0 . 



GC TO TOO 
00 l F ( W F . L T . 3 2CC . ) GO TO 3 



TFTwFTGTTTZJOUTT 

lfj.=pf 

..GC; TP ?QQ 



GU 10 4' 



WF1. =3300 . 
GC TO 200 
■■ wri-l?.300. 



.00 i = l F T X( X3/1C3J. J+l 

J= T F T X ( (WF 1-3300. ) /TOQO. ) + l 



~C tw = V "3 PITT 

DW=’aFI-VvFT ( J ) 

-O ELO LU a nR /l O OO.) » tC ET ( :- U «-J 4 --=flET (-U-LU. 



OEl fjw*(PW/lC00. )*(GET(I,J+L)-GET( 1 * J ) ) 
T A6QE=QET { I , J J+OEL CR+OEL OW 



■■ k r T l if< » . ; 



exn 



T>‘F FOLLOW IPG OATA ARE PART OF THE PLOTTING P A C F A G r 



FPk USL/360 S i MUL A 1 10 
SHOIJ! 0 BF MCfJTFTEn T F 

— LS M LA lLAhl F 



. LANGUAGE A I M P S . — FFTpr 
OTHER THAN THIS PLOTTING PACKAGF 
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'Plot .5Y r .T M r j 1) * 


;FI.Ai\r*NT c 1. fi 

J] L* =-2_. 2_. 




C 




tr Cf- T HP- V ART AT TOM 0 F PROPELLER SPEFC 
>LA.\T:NfS IA 

ps — a a o , } c c , 


VS .TIME 

z 


c; 




.CT CF °pfil , PPO'l , PERCl , PFRC2 ,P WAX 
, FL AN T J N T S 1.A 




r*-* 




nj u . i h u . U 6 u . 

GT CF 7 A ANO VS.TT^F 

'.f\ ;lM~P r S 1 A 


5T 


T: 




,U 1. 0.0 150. 

l T CF R8 VS.TIVE 

; r L >T T • A ! 1 S 1-4 — 


5 . 


5. 


• 



,o' ' 0.38 0.0 0.004 5 . 5. 
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